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Abstract. A successful description of hadron-hadron collision data demands a profound un¬ 
derstanding of quantum chromodynamics. Inevitably, the complexity of strong-interaction 
phenomena requires the use of a large variety of theoretical techniques—from perturbative 
cross-section calculations up to the modelling of exclusive hadronic final states. Together 
with the unprecedented precision of the data provided by the experiments in the first run¬ 
ning period of the LHC, a solid foundation of hadron-hadron collision physics at the TeV 
scale could be established that allowed the discovery of the Higgs boson and that is vital for 
estimating the background in searches for new phenomena. This chapter on studies of quan¬ 
tum chromodynamics at the LHC is part of a recent book on the results of LHC Run 1 [Tj 
and presents the advances in theoretical methods side-by-side with related key measurements 
in an integrated approach. 
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1 Introduction 

Quantum chromodynamics (QCD) is the well-established quantum field theory of the strong interaction [2, 
g] and one cornerstone of the Standard Model (SM) of particle physics. Like the electromagnetic and the 
weak force, QCD belongs to the class of gauge field theories. The assumption that the corresponding 
gauge symmetry is an exact symmetry of nature results in truly massless force carriers of the strong 
force, the gluons. However, as a consequence of the non-Abelian character of the SU(3) QCD gauge 
group, the gluons carry a strong or so-called “colour” charge and interact amongst themselves. This is 
a striking difference to the electromagnetic force, which is mediated by electrically uncharged photons, 
and it induces—amongst other effects—the confinement of strongly interacting particles at low energies. 
In the following, the colour-charged constituents of hadrons, i.e. quarks, antiquarks, and gluons, will 
generically be denoted as “partons”. 

In hadronic collisions like at the LHC, QCD effects are omnipresent, and their detailed understanding 
is indispensable for the interpretation of collider data, whether to search for new phenomena or to perform 
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precision studies of model parameters. Despite complications through the gluon self-interactions, the per¬ 
turbative expansion of QCD (pQCD) that describes interactions with large momentum exchanges (“hard 
interactions”) in terms of parton-parton scatterings remains the most powerful theoretical technique. This 
technique allows quantitative predictions at parton level for observables ranging from inclusive produc¬ 
tion rates to shape observables, which are sensitive to details of the QCD radiation pattern. Nevertheless 
the theory’s confining nature obliges us to address non-perturbative aspects. This includes a reliable 
understanding of the short-distance parton structure of the initial-state protons, the fragmentation of 
final-state partons into hadrons, or the modelling of soft proton interactions. 

The unprecedented experimental precision achieved with the new detectors at the LHC requires 
equally accurate theoretical predictions and has sparked rapid progress in the field of perturbative cal¬ 
culations using both analytical methods and modern Monte Carlo (MC) event generators. The interplay 
between experiment and theory enforced the development of new observables and novel techniques to 
match the challenges arising on both sides. In summary, the improved theoretical understanding of the 
QCD dynamics including the ability to precisely predict even complicated high-multiplicity final states 
and the excellent performance and understanding of the LHC machine and detectors together with so¬ 
phisticated analysis techniques reveal a more refined and detailed picture of QCD than ever before. 

For the presentation of QCD-related experimental results obtained during LHC Run 1 and the under¬ 
lying theoretical developments, an integrative approach is chosen—reflecting the productive and fruitful 
interplay of the two communities. As a consequence, compromises on the content had to be taken, and only 
a selection of the most important measurements is discussed—omissions in the presentation of theoretical 
methods were unavoidable. 

The chapter at hands is structured as follows: After a brief reminder of the basics of the QCD theory 
and the central aspects of perturbative QCD, the discussion focuses on various approximations in the 
modelling of scattering processes, including parton-shower simulations and parton-level predictions at 
next-to-leading and next-to-next-to-leading order in the strong coupling. Then the discussion turns to 
the observation of nmlti-jet final states, the successful description of which requires the combination of 
both fixed-order calculations and all-order (i.e. parton-shower) techniques. This is followed by a section 
on analytical methods for the resummation of large logarithms; these are exemplified using gap-fraction 
and jet-substructure observables. The chapter closes with a presentation of various phenomena and mea¬ 
surements sensitive to non-perturbative aspects of QCD. 


2 Basic Elements of QCD 


The discussion shall begin with a brief reminder of the ingredients of the QCD Lagrangian that defines the 
Feynman rules required for a perturbative analysis of QCD. The classical QCD Lagrangian is composed 
out of the free Dirac Lagrangians for the six quark fields and the kinetic and self-interaction terms for the 
gluon fields A“, labelled by a colour index a = 1,..., 8. These two parts get minimally coupled through 
a gauge covariant derivative 

£qcd = £ gauge H" -^quarks ? 


where 


r — _ i p a pw 

■‘■'gauge — | r p r o : 


with = d,Al - d v Al - g s fabcA b A c u 


the gluon field-strength tensor, and 


£quarks = 9 (*7 ^ D U ~ ™q) Q , with D^ + ig s t a A“ 

q£{u,d,s,c,b,t} 


the QCD covariant derivative (see also the introduction to the SM Lagrangian in Chap. 4 of Ref. CD- 
Quark masses are denoted by m q . The SU(3) generator matrices introduced here obey the algebra 

[t a ,t b ] = ifabct° , 

defining the QCD structure constants f a b c • The classical QCD Lagrangian exhibits the property of local 
gauge invariance, i.e. invariance under a simultaneous redefinition of the quark and gluon fields. As 
a consequence of this internal symmetry, it is impossible to define the gluon field propagator without 
explicitly specifying a choice of gauge. A Lorentz-covariant way to fix the gauge is given by the class of 
R,£ gauges, imposed by adding a term 
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to the classical Lagrangian. Because of the non-Abelian character of the QCD gauge group, the full 
Lagrangian of the quantum field theory features a further contribution, the ghost Lagrangian 

£ghost = d^rj t (D^ b ri ) , 

that represents the field-dependent Faddeev- Popov determinant. The ghost fields r] a are represented 
by anti-commuting scalar fields. This completes the Lagrangian for a consistent version of a quantum 
field theory of the strong interaction. Accordingly one can read off the QCD Feynman rules, featuring 
three-point quark-quark-gluon and ghost-ghost-gluon interactions as well as triple and quartic gluon 
self-interactions. All of these interaction vertices are proportional to the strong charge g s . This is also 
the relevant parameter when applying the method of perturbation theory to QCD. Defining the QCD 
counterpart of the QED fine-structure constant a s = g 2 /4n, one can expect a truncation of the power- 
series expansion for a given observable O , i.e. 

O = Of) + (D\(X S + 0 2 OC s + ... , 


to yield meaningful estimates as long as a s -C 1. 

A prime example of a quantity evaluated in perturbation theory is the QCD j3 function. It determines 
the running of the coupling constant a s through the renormalisation group equation 


Q- 


da s 
dQ 2 


= (3(a s ), with /3(a s ) =-a 2 s (b 0 + bxa s + b 2 a 2 + (D(a 3 s )), 


and 


bo = 


33 — 2 tif 


bi = 


153 — 19tif 


bo = 


77139 - 15099n F + 325n|, 


(1) 


( 2 ) 


12?r 7 24t r 2 ’ 3456 tt 3 

tif denotes the number of quark flavours with masses m q smaller than the scale Q. Note that the 
higher coefficients b 2 and 63 (see Ref. m are renormalisation-scheme dependent. Here b 2 is quoted in the 
MS scheme. Retaining only the leading term bo, equation ([TJ is solved by 


a s (Q 2 ) 


_ Usjg 2 ) _ 

1 + 6 0 In (Q 2 / g 2 ) a s (/j, 2 ) 7 


( 3 ) 


which relates the strength of the coupling at a scale Q to the one at scale g, assuming both scales to be 
in the perturbative regime. The non-Abelian nature of QCD manifests itself in the negative sign of the (3 
function. Thus, as long as uf < 17, the coupling becomes weaker at higher scales Q, or, in other words, the 
QCD colour charge decreases when the distance decreases. For high scales Q, QCD becomes almost a free 
theory—a property known as “asymptotic freedom”. It is this weakly coupled regime where perturbative 
methods can successfully be applied and quantitative predictions for hard scattering processes can be 
made. The world average value of the strong coupling as of 2014, quoted at the scale of the Z-boson mass 
Mz, is given by 

a s (M z ) = 0.1185 + 0.0006, 

derived from hadronic r-lepton decays, lattice QCD calculations, deep-inelastic scattering data, electron- 
positron annihilation processes, and electroweak precision fits [5]. Figure [l] shows a summary of measure¬ 
ments of the strong coupling at energy scales ranging from the mass of the r-lepton of M T « 1.8 GeV up 
to the TeV scale thanks to newly included LHC data. The historical development of a s determinations 
is discussed in Chap. 12, Fig. 12.2, of Ref. [Tj. 

The dynamical behaviour of a s (Q 2 ) implies an increase of the QCD coupling at small momentum 
transfer, i.e. large distances. When the coupling approaches unity, perturbation theory is not valid any¬ 
more. The parameter Aqcd is defined as the scale, where a s (Q 2 ) formally diverges. With this definition, 

equation ([3]) can be rewritten as a s (Q 2 ) = (6 0 ln (Q 2 /A qc D )) 1 . For nj? = 5 flavours in the MS scheme, 
Aqcd roughly amounts to 214 MeV and represents the dividing line below which one is in the manifestly 
non-perturbative regime of QCD. It is the growth of the coupling at small scales that makes QCD a the¬ 
ory of the strong interaction—the fundamental force that confines the quarks and gluons into ordinary 
hadronic matter, e.g. the protons and neutrons. For the purpose of LHC physics, one has to account for 
this phenomenon of “confinement” when modelling the transition from free quarks and gluons to the 
bound-state hadrons observed in the detectors. Lacking a first-principles understanding of this process, 
one mostly has to rely on MC models for this aspect. Even more fundamentally, the partonic content of 
the colliding protons needs to be parametrised in order to allow for a description of LHC collision events 
through partonic scattering processes. 
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Fig. 1. Summary of measurements of the strong coupling a s as a function of the respective energy scale Q. By 
including new LHC data, the range in Q could be extended to the TeV scale. (Adapted from. Ref. f "5j/. ) 


3 Perturbative QCD 

Quantitative predictions based on the non-Abelian QCD Lagrangian can be obtained either with lattice 
methods employing a discretised space-time or using perturbation theory. Given the complexity of the 
final states produced in proton-proton collisions with high momentum transfer, lattice techniques are of no 
practical importance for the prediction of LHC events at present. Instead one has to rely on perturbative 
methods, which are possibly supplemented by models for the transition of partons to hadrons. 


3.1 Cross-Section Predictions 


In perturbative QCD, the cross section for a hard scattering process at a hadron-hadron collider can be 
written in the following factorised form 


&h 1 h 2 ^x — 



dx±dx2 fi/H 1 {xi,HF)fj/H 2 {x2,HF) 


X&ij-tX 


(x 1 P 1 ,x 2 P2,ci s (rr), 



( 4 ) 


where the sum extends over all contributing initial-state partons i,j G {q,q,g}. It is assumed here that 
the scale associated with the hard process, Q , is much larger than Aqcdj the delimiting scale for the 
applicability of perturbative methods to QCD. In this high-energy limit, effects related to the binding 
of the partons in the initial-state protons can be neglected. As a consequence, the cross section simply 
factorises into a product of parton distribution functions (PDFs) Rf) of non-perturbative origin 

and the perturbatively calculable partonic cross section. The parton distribution functions model the 
probability to find a parton of flavour i in the incoming hadron H (protons at the LHC) with a fraction x 
of the hadron’s momentum P. The resulting squared partonic centre-of-mass energy is given by s = X\X 2 s, 
with s = (Pi + P 2) 2 the squared hadronic centre-of-mass energy. 

The collinear factorisation ansatz underlying equation (El) is the key to quantitative predictions in 
the framework of QCD that can be compared to actual LHC collision data. Based on the property of 
asymptotic freedom of QCD, the desired cross section can be expanded as a power series of the coupling 
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constant a s . The lowest-order coefficient is denoted as leading order (LO), the subsequent ones as next- 
to-leading order (NLO) and next-to-next-to-leading order (NNLO), respectively. 

When calculating the LO, NLO or NNLO estimate for a partonic cross section, the appropriate QCD 
evolution of the parton distribution functions has to be used. 

Leading-order cross-section calculations are fully automated by means of tree-level matrix-element 
generator programs such as Alpgen [B], Amegic [7j, Comix [5], Helac/Phegas [5J, MadGraph m or 
Whizard m - These codes are capable of providing integrated cross sections and parton-level events for 
almost arbitrary Standard Model final states, with multiplicities ranging up to ten particles. In particular 
for high-multiplicity final states, implementations relying on recursive algorithms for the generation of 
the expressions for the amplitudes, e.g. Berends-Giele recursion m, prove most efficient [Ullll- 

Over the past years there has been enormous progress in the evaluation of processes at NLO and 
NNLO in the strong coupling. These developments and related precision measurements will be addressed 
in detail in later sections of this chapter. 


3.2 Fragmentation and Hadronic Jets 

One entity from equation p} that has not yet been discussed is the final state A of a collision. The simplest 
reaction that can be considered is the Drell-Yan process m, where a quark and an antiquark annihilate 
to produce a lepton pair: a ( qq —► i + t ~). In this case there are no strongly interacting particles in the final 
state, and the theory prediction can directly be compared to the measured leptons. Merely the proton 
remnants, which fragment into hadrons along the beam lines, have to be described by non-perturbative 
models in MC event generators. At high transverse momenta, the two leptons are well separated from 
any such proton debris and high-precision comparisons with theory even at NNLO become possible. This 
is discussed in more detail later in this chapter. 

However, in the vast majority of reactions at least some colour-charged partons are produced so that 
a further step covering the transition from the partonic final state to measurable particles, the so-called 
“particle level”, is needed. Here, “measurable” refers to colour-neutral particles with mean decay lengths 
such that ct > 10 mm, where c is the speed of light and r the lifetime of a particle. One possibility to 
account for this transition is to reuse the concept underlying the PDFs that describe the partonic content 
of a hadron, only in an inverted sense. The necessary functions Dk->h(z, hf) are called fragmentation 
functions (FFs) and are the final-state analogues of the PDFs. They parametrise the probability of 
finding a hadron h within the fragmentation products of parton k, carrying the fraction z of the parton 
momentum. Like the PDFs, fragmentation functions depend on a non-physical resolution or fragmentation 
scale fj,p- Again, these functions can currently not be determined by first principles in QCD, but once 
they have been measured (for example under the experimentally more favourable conditions of e + e - 
collisions at the LEP collider), they are universally valid. 

A second possibility to account for the transition to measurable particles makes use of the concept of 
hadronic jets. Instead of looking into the detailed production of identified particles—an experimentally 
very challenging endeavour- for the majority of processes it is sufficient to know how much energy and 
momentum is carried away by hadrons. QCD predicts that large-distance non-perturbative (NP) effects 
are mostly decoupled from the hard reaction so that highly energetic partons fragment into a collimated 
stream or “jet” of hadrons, which inherits energy and momentum from its parent parton. To define what 
“collimated” means, a prescription is required that, given some distance measure, unambiguously decides 
which objects belong to a jet. As one wants to compare predictions by pQCD with measured particles, 
tracks, or energy depositions as illustrated in figure [2j a jet algorithm is needed that is applicable to 
theoretical calculations as well as to measurements from different experiments. 

Most importantly, to deal with the cancellation of collinear and soft singularities appearing in pQCD, 
a jet algorithm must be collinear- and infrared-safe. This means that the outcome of a jet-clustering 
procedure depends neither on the splitting or merging of collinear parton four-vectors nor on the addition 
of arbitrarily soft partons to the list of objects to be clustered. The first description of a collinear- and 
infrared-safe jet algorithm, which grouped partons or particles together that are inside an angular cone 
around a specific direction, was given by G. Sterman and S. Weinberg in 1977 HU- 

The most important requirements for jet algorithms from the experimental side are i) independence of 
detector details, ii) maximal reconstruction efficiency, iii) minimal resolution smearing, iv) computational 
efficiency, and v) ease of calibration. The extension of the original jet definition, which was specialised to 
2 -jet events in e + e _ collisions, to all kind of reactions and the partially conflicting requirements lead to 
various new propositions that were tried and discussed in several workshops HZHI2|. Two classes of jet 
algorithms emerged: 

1 . cone algorithms that geometrically assign objects to the leading energy-flow objects in an event; 
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Fig. 2. Illustration of a jet to which bundles of partons, hadrons, or detector measurements are grouped together. 


2 . sequential-recombination algorithms that iteratively combine the closest pairs of objects [20] . 


Many cone algorithms need starting points with a minimum energy or momentum for the cone 
directions—so-called “seeds”. These spoil the condition of collinear safety. A seedless infrared-safe cone 
algorithm exists in the form of the SisCone algorithm {21], which avoids this problem. For reasons of 
computational efficiency, though, the method of choice employed at the LHC is the anti-fc t clustering 
algorithm [22] as implemented in the Fast Jet package [23]. For the jet size parameter R (the equivalent 
to the cone size in cone algorithms) the LHC collaborations chose the values 0.4 and 0.6 for ATLAS 
respectively 0.5 and 0.7 for CMS. An extensive overview of jet definitions in QCD and their history is 
presented in Ref. [231 . 

The most fundamental quantity of jet production that can be investigated is the inclusive jet cross 
section as a function of the jet px and rapidity y, d 2 a/dp!^ dy, where every jet of an event contributes. 
Figure [3] presents a summary of inclusive jet cross section measurements, performed with various jet 
algorithms, from pp or pp collisions at centre-of-mass energies from 546 GeV up to 7 TeV at central 
rapidity. The measured cross sections stretch over eleven orders of magnitude and range up to 200 GeV 
at SPS, 600 GeV at the Tevatron, and 2 TeV at the LHC. 

LO predictions (not shown) correctly describe the steep drop with increasing p^ 4 , but fail to predict 
the absolute values more precisely than in an order-of-magnitude estimation. More accurate calculations 
at NLO are able to improve this situation. 
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Fig. 3. Summary of inclusive jet cross-section measurements at pp or pp colliders. Data taken from Refs . [25L[35l 


4 Parton Showers: The Bulk of the Emissions 


The strength of the collinear factorisation ansatz allowing us to define universal PDFs and FFs lies in its 
repeated applicability. It enables us to define evolution equations such as the DGLAP equations , 

which determine the variation of the PDFs and FFs when the factorisation scale /ip is changed. For the 
case of the initial state PDFs one derives 


2 9fi(x 


Af) 


d/4 


E 

i={9,9.9} ; 




(5) 


with Pij the (plus-prescription) regularised Altarclli Parisi splitting functions 

= c r (DAzili) , p„w = C F (d±A + ? {(1 
P qg {z)=T R (z 2 + {l-z 2 )) , 

p gg( z ) = 2Ca H — F z{ 1 - 2 )^ + <5(1 - z) 



- 4 n F T R 
6 


( 6 ) 


Furthermore, it should be noted that P gq = P gq and P qq = P qq . While the LO, i.e. one-loop, approximation 
is quoted here, the QCD splitting functions are known up-to NNLO, i.e. three-loop accuracy HUE]. 
The LO DGLAP evolution allows for a interpretation by means of simple branching processes. A parton 
i resolved at scale p, F may have originated from a branching of parton j resolved at some higher scale. 
This transition of parton j to i is accompanied by the emission of an additional QCD parton. When 
applying the DGLAP equations to solve for the scale evolution of PDFs or FFs, these emitted particles 
get ignored, by considering inclusive processes only. 

However, in parton-shower MC programs these emissions are made explicit, and the subsequent 
branching of initial-state and final-state partons results in a cascade-like picture, modelling the initial- 
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Fig. 4. Measured normalised /3 distribution, corrected for detector effects, in comparison to various MC pre¬ 
dictions. The error bars indicate the statistical uncertainties, while the shaded bands correspond to the total 
systematic uncertainty. (Adapted from Ref. ]SOj.) 


and final-state evolution of the hard process’ parton configuration into an ensemble of many QCD par- 
tons. In essence these parton-shower emissions approximate higher-order real-emission corrections to the 
hard scattering process. Shower simulations form an integral part of MC event generators—simulation 
tools to describe the fully exclusive hadronic final states of individual scattering events [52]. They link 
the hard scattering at some high scale Q 2 to an exclusive set of partons with typical separation scales 
of order Qq ss 1 GeV 2 , the cut-off scale of the shower evolution. Through this parton-shower link it is 
possible to invoke universal models for the hadronisation of partons into hadrons, independent of the hard 
process scale Q 2 . Physically, parton showers account for the intra-jet evolution of jets by accounting for 
the emission of almost collinear partons, and they model the inter-jet QCD activity through the emission 
of soft, wide-angle partons. 


4.1 Colour Coherence 

In particular the correct incorporation of soft-gluon emissions needs special consideration. In contrast to 
collinear emissions, which factorise at the cross-section level, soft-gluon emissions factorise at the level of 
individual QCD amplitudes. One should therefore consider soft gluons to be emitted by the scattering 
process as a whole, given by the squared sum of all contributing amplitudes. At first glance this approach 
seems to spoil the parton-shower picture of independent subsequent emissions. However, soft colour- 
coherence effects originating from the interference of individual amplitudes can be incorporated into 
parton-shower simulations by suitable choices for the shower-evolution variable. Most prominently this 
can be achieved by using the emission’s opening angle as ordering parameter, as employed in the Herwig 
and Herwig++ generators [431 144]. 

However, angular ordering is not the only option to include colour-coherence effects. A lot of effort 
went into the construction of new parton-shower algorithms based on subtraction formalisms as they 
are used in NLO QCD calculations. These new shower schemes, either based on Catani-Seymour dipole 
factorisation pToHTTl or antenna subtraction 0SJ , implement soft-gluon coherence through a partial 
fractioning of the QCD radiator antennas, thereby introducing the notion of an emitter, i.e. splitter, and 
an associated spectator parton that accompanies the splitting process. 

In order to study QCD coherence effects one needs to devise observables sensitive to rather soft 
emissions. This can be achieved by selecting final states with at least three jets exhibiting a sizeable 
spread in transverse momentum between the hardest and the softest jet. One such analysis was presented 
by the CMS collaboration in Ref. m, based on a data set with an integrated luminosity of 36 pb 1 
collected in 2010. The analysis inspects events with at least three anti-A;* jets using a distance parameter 
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of R = 0.5, ordered in transverse momentum such that px,i > px ,2 > Pt, 3 • The event selection criteria 
are given by 


Pt,i > 100 GeV, pt ,3 > 30 GeV, |t 7 i i2 | < 2.5 
Mi 2 > 220 GeV, 0.5 < AR 23 < 1.5. 

The observable considered to probe colour-coherence effects, called /3, is defined as the azimuthal angle 
of the third jet with respect to the second jet, i.e. 

tan/3 = „ 

At/23 

In the presence of QCD coherence, the emission of the parton initiating the third jet is expected to 
preferentially he in the event plane defined by the emitting parton and the beam axis. Figure [4] shows 
the normalised /3 distributions measured for two regions of the second-jet pseudo-rapidity, i.e. central 
1 772 1 < 0.8 and forward 0.8 < I772 1 < 2.5. The observable is thereby defined as -F?7 2 ,i(/3) = iV^/A^, with N v 
the total number of events in the respective 772 region and N rh . c the number of events in the ith /3 bin. The 
data are compared to various particle-level MC predictions, i.e. specific tunes of Pythia6, Pythia8 m 
and Herwig 4—b, all based on LO 2 —)• 2 matrix elements plus parton showers, and a combination of 
MadGraph using exact 2 — > 2 and 2 — > 3 tree-level matrix elements supplemented by PythiaO parton 
showers. It has been shown in Ref. m that the data clearly support the inclusion of coherence effects 
in the simulation; the regions of small /3 and /3 ss 7 t/2 are particularly sensitive. However, none of the 
generators used in the analysis describes the data satisfactorily over the entire phase space. Herwig-I —f- 
models the data best but the agreement is rather poor in the forward region around (3 ~ 7r. The inclusion 
of the exact 2 —> 3 matrix element from MadGraph slightly improves the pure Pythia6 parton-shower 
prediction. 


4.2 Azimuthal Decorrelation 

Clearly, the inclusion of parton showers is indispensable when attempting to simulate event configura¬ 
tions sensitive to the emission of multiple partons. Exemplary quantities are event shapes |52l454| . or 
sub- or intra-jet-related observables [55H58] . Here, the measurement of dijet azimuthal decorrelations as 
performed by ATLAS and CMS [59j :60] shall briefly be discussed. The analyses evaluate 36 pb -1 and 
2.9 pb -1 of pp collision data at yfs = 7 TeV by ATLAS and CMS, respectively. The observable considered 
is the azimuthal angle Acf> between the two jets leading in jet px- For events with exactly two high-px 
jets and nothing else, the azimuthal angle between the two jets is fully correlated through momentum 
conservation in the transverse plane and Acf) ss it. Multi-jet production disturbs this balance so that 
the two leading jets become decorrelated in azimuthal angle and A</> < 7 r. Hence, the observable Acf) is 
indicative of multi-jet production, while at the same time only the azimuthal angles of the two leading 
jets are measured, avoiding the large uncertainties associated to the jet energy calibration. 

With Acf> it is possible to probe complementary aspects of perturbative QCD. The region A(j> ss 7r 
is sensitive to multiple soft and collinear emissions and thus requires a proper resummation of soft 
and collinear logarithms. The region of large azimuthal decorrelations, i.e. Acj> <C 7r, indirectly probes 
the production mechanism for additional hard jets and requires the inclusion of higher-order matrix 
elements. Figure [5j a) shows a comparison of the normalised distribution of Acf) for five slices in leading- 
jet px as measured by CMS [BO] to perturbative predictions for three-parton production at LO and NLO, 
respectively, in the region 2tt/ 3 < Acf) < 7r. The leading order has predictive power only in a very narrow 
range of Acj). The NLO calculation improves this substantially, but fails when approaching Acf) ss 7r, and 
below 27 t/ 3 matrix elements for 4-jet production are required. 

In figure [5jb) the normalised Acj> distributions as measured by ATLAS |59j are compared for nine 
slices in leading-jet px to predictions from the MC event generators Pythia6, Herwig, and Sherpa. 
Pythia6 and Herwig are based on LO 2 —> 2 QCD matrix elements matched with parton showers. 
Sherpa additionally includes tree-level matrix elements for 2 —> 3-6 jet production properly matched to 
its parton-shower algorithm mm- All three generators describe the data well over the measured range 
of 7r/2 < Acf> < 7 r and all p“ ax slices. This includes both the region ir/2 < Acf> < 2ir/3 where multi-jet 
contributions are significant and Acf) ss 7 r that cannot be described by a fixed-order calculation. 
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Fig. 5. Normalised differential cross section (1 //AA(j>) for multiple slices in leading-jet px: (a) CMS data 
compared to fixed-order predictions in pQCD; (b) ATLAS data compared to predictions by the MC event gener¬ 
ators Pythia6, Herwig, and Sherpa. (Adapted from auxiliary material provided with Refs. [5ffl [6$ .) 


5 NLO: The New Standard 


For most observables, NLO predictions in the strong coupling represent the first accurate theoretical 
estimate that allows an assessment of associated theoretical uncertainties. The evaluation of NLO cross 
sections, however, is more involved and requires to consider real-emission and virtual one-loop corrections 
that are individually singular in the infrared region. While these divergences cancel in the sum of both 
contributions for sufficiently inclusive observables, they render the numerical evaluation of such cross 
sections difficult. Several state-of-the-art techniques exist that allow the exact cancellation of the diver¬ 
gences separately for the real and virtual contributions through the introduction of suitable subtraction 
terms. Most widely used is the dipole-factorisation method by Catani and Seymour [H2J; alternatives are 
provided by the Frixione-Kunszt-Signer (FKS) method [B3|, antenna subtraction [B2], or the recently 
developed Nagy-Soper formalism [65 : , f66]. 

The enormous progress recently experienced in the field of NLO QCD calculations was sparked by two 
important developments: i) the introduction of fast and efficient methods for the calculation of virtual 
amplitudes, see for instance Refs. EM2; ii) the organisation and implementation of complete NLO 
calculations in the framework of parton-level Monte Carlo event generators such as Helac/Phegas [73] . 
MadGraph5_aMC@NLO [73], or Sherpa [75a 76]. All these approaches employ automated subtraction- 
term generators, see Refs. E51I33MII, that implement the Catani-Seymour or FKS subtraction formalism. 
The real-emission corrections as well as the phase-space integration are handled by tree-level matrix- 
element generators such as Amegic [7], Comix [3], MadGraph [TU], MadGraph5_aMC@NLO [73] 
EU, or Helac jJJ]. Virtual amplitudes, typically provided by specialised one-loop generators such as 
BlackHat [72], Gosam [82], Helac-ILoop [83], MadLoop [84], Njet [85], OpenLoops+Collier [86J 
or Recola m can be incorporated via the universal BLHA interface [55J [39] or dedicated solutions. 
Examples of recent NLO calculations that have been performed using a combination of the tools listed 
above include: W + 4,5 jets [Ml 01], Z + 4 jets [92] , 4-jet and 5-jet production [93H55] , tt + 2 jets [SB] and 
77 + 3 jets m- Most of these new tools are readily available to perform NLO QCD event generation for 
use in LHC data analyses. 
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5.1 Jet Counting 

With NLO predictions at hand, more precise comparisons of data and theory become possible. The 
production of single inclusive jets and of dijets constitute the most basic QCD processes at hadron 
colliders. Figure [g] shows a comparison between data and theory at NLO, i.e. 0{c\a 2 s + C 2 a 3 ), for (a) 
inclusive jets from CMS 01] and (b) dijets from ATLAS 05], both at 7 TeV centre-of-mass energy. 
Recently, also trijet cross sections were measured E3 nmj. The inclusive jet cross section has been 
measured as a function of and y in ranges of 0.1 < px < 2.0 TeV and \y\ < 2.5 (see Refs. 03 i33])• 
The dijet cross section is presented as a function of the dijet mass, TO 12 , and rapidity separation, y* = 
\yi — 2 / 2 1 /2, of the two highest-px jets in regions of 0.3 < TO 12 < 5.0 TeV and y* < 3.0, respectively. 
In both cases the NLO theory has been corrected for non-perturbative effects, and for the more recent 
ATLAS publication also electroweak (EW) corrections have been considered. While EW effects start being 
sizeable at high p^ and \y\ or large mi 2 lOl . non-perturbative effects, estimated via their modelling 
in MC event generators, become negligible. The data are well described by theory over many orders of 
magnitude in cross section and in wide kinematic ranges so that fits of parameters entering the pQCD 
calculations become feasible. 
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Fig. 6. (a) Inclusive jet cross section from CMS and (b) dijet cross section from ATLAS, both at 7 TeV centre- 
of-mass energy, in comparison to theory at NLO including non-perturbative (CMS) or non-perturbative and 
electroweak corrections (ATLAS). The anti-fct algorithm has been applied with the jet-size parameter R being 
0.7 for CMS and 0.6 for ATLAS. (Adapted from Refs. I9H \98 ] i .) 


5.2 Jets and the Gluon PDF 

In the following, PDF fits are considered as a first example of determining parameters entering pQCD 
calculations. Many choices are possible in terms of data set selection, theoretical ingredients and approxi¬ 
mations, parameterisation, or fitting techniques and criteria. There is no unique approach to PDF fits. The 
PDF sets in use at the LHC are produced by six PDF fitting groups: ABM, CTEQ, (G)JR, HERAPDF, 
MSTW, and NNPDF (10214107] . where the quoted publications refer to pre-LHC determinations. The 
fundamental experimental input central to all PDF determinations are deep-inelastic scattering (DIS) 
data collected by the HI and ZEUS collaborations at the ep collider HERA. The HERA collaborations’ 
own PDF fits based solely on their data are presented as HERAPDF sets. The ABM group extends 
the choice of data to other DIS measurements and fixed-target Drell Yan data, while the (G)JR group 
also considers high-px jet production, which cannot yet be described at NNLO, to better fix the gluon 
PDF. The trio of CTEQ, MSTW, and NNPDF finally try to incorporate as many data sets as possible 
and to provide truly global fits. New PDF sets including LHC data became available recently or are in 
preparation for the next start-up of the LHC with 13 TeV center-of-mass energy. 
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It is possible to judge the potential impact of new data, e.g. from the LHC, on a given PDF set without 
a complete refit. The technique of “Bayesian reweighting” was first implemented in the context of the 
NNPDF approach 1108 !. The NNPDF collaboration |i l05j provides PDF sets in the form of an ensemble of 
replicas, which sample variations in the PDF parameter space within allowed uncertainties. Their method 
fits a PDF to smeared data points, where each data point is sampled from a multi-Gaussian distribution. 
The central value of the distribution is equal to the measurement point and the covariance is taken from 
the experimental uncertainties. This procedure is repeated N times and provides an ensemble of N PDFs 
representing the uncertainty of the PDF fit. Hence, the NNPDF prediction for an observable is given by 
the mean and standard deviation as estimators for the true value and its uncertainty. Via the Bayesian 
reweighting, the impact of new data can be approximately taken into account by giving according weights 
different from unity to each replica. If this leads to many replicas with zero or very small weights, this 
technique cannot be applied and refits must be performed. Otherwise the re-evaluation of an observable 
while using the re-weighted set of replicas approximately accounts for changes in the PDFs due to the 
additional data. 

In a similar way to the mean and standard deviation as estimators for the main prediction and 
uncertainty of an observable, the NNPDF ensemble of replicas can be used to investigate the correlation 
between an observable and the PDFs as a function of energy scale Q and momentum fraction x. As an 
example, figure[7[a) shows the correlation coefficient p between the inclusive jet cross section as measured 
by CMS m and the gluon PDF. For x > 0.01 a large positive correlation is exhibited between the two. 

In case of the other five PDF sets, the method of diagonalisation of the Hessian matrix [ 111] (also 
called “eigenvector method”) is employed to express uncertainty estimates for their respective choices of 
PDF parameterisation. Their prediction for an observable in pQCD can thus be evaluated from one central 
PDF set with an uncertainty given by quadratic addition of all deviations obtained while using additional 
PDFs that correspond to variations along the directions of the eigenvectors. A Bayesian reweighting was 
shown to be feasible here as well m, and an alternative approach is studied in Ref. m- 
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Fig. 7. (a) Correlation coefficient p between the inclusive jet cross section and the gluon PDF as a function 
of energy scale Q and momentum fraction x. (b) PDFs as determined from fits to HERA DIS data alone or in 
combination with CMS jet data. The uncertainties on the gluon PDF at high x are reduced compared to a fit 
with DIS data alone (uncertainties not shown). (Adapted from Ref. j 110/ and auxiliary material.) 


Employing an open-source tool like the HeraFitter project 105 . ! 114l lll5i . which provides fitting 
code, data sets, and theory predictions in a common framework, it is even possible to perform complete 
PDF fits. The CMS collaboration used this tool to study the impact of their inclusive jet measurements 
on the gluon PDF [110] following the fitting setup and using the same DIS data as described in Ref. 11051 . 
The shape of the PDFs at a starting scale Qq is assumed to be: 
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xg(x) = A g x Bg (\ — x) c9 — A' g x B »(1 — x) c s , 
xu v (x) = A Uv x B u ^ (1 - x) Cu * (1 + E Uqj x 2 ), 
xd v {x) = Ad„x Bd '' (1 — x) Cdv , 
xU(x) = Ajjx B ij (1 — x) c u , 
xD(x) = Ajjx b °(1 — x) c ^. 

Here, xg(x) is the gluon distribution, xu v {x ) and xd v (x) represent the valence quarks with u v (x) = 
u(x) — u(x), d v (x) = d(x) — d(x), and U(x), D(x) are the up- and down-type antiquark distributions 
U(x) = u(x), D(x) = d{x) + s(x). The indexed symbols Ai, Bi , Ct with i £ { g,u v ,d v ,U,D }, A' , B', 
C' g , and E Uv stand for 19 parameters in the definition of the PDFs. Note that there are no heavy quarks 
at a starting scale that is chosen to be Qq = 1.9 GeV 2 . The evolution of PDFs in Q 2 follows from QCD 
according to equation ([5|), where in the HeraFitter setup a generalised-mass variable flavour number 

scheme uminzi is used. 

A g , A Ufu , and A ^ are normalisation parameters that are constrained by QCD sum rules. The B and 
C parameters describe the limiting behaviour when very small, x —> 0, or very large, x —¥ 1, proton 
momentum fractions are approached. Additional terms with parameters for the gluon and it-valence 
distribution allow for some more flexibility in shape, while constraints of Bjj = Bj^ and Ajj = Aq(1 — f s ) 
are applied to ensure the same normalisation for the u and d densities at high x. The strangeness fraction 
is set to f s = 0.31, as obtained from neutrino-induced dimuon production 1 118] . The parameter C’ g is fixed 
to 25 [ 10611117] and the strong coupling constant to a s (Mz ) = 0.1176 so that in total 13 free parameters 
remain to be determined. 

The PDF fit result using HERA DIS data alone or in combination with CMS jet data is presented 
in figure [7j]b). As expected, a clear impact on the gluon PDF for momentum fractions x > 0.01 and a 
reduction of uncertainties (not shown) is exhibited. 
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Fig. 8. The ratio of the inclusive jet cross sections at yfs = 2.76 and 7 TeV as measured by ATLAS divided by 
a pQCD prediction at NLO (with CT10 PDFs) of the same ratio. Experimental uncertainties are indicated by 
error bars (statistical) and gray boxes (systematic). The uncertainties inherent to the PDFs are represented by 
bands. (Adapted from Ref. I35\l . ) 
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5.3 Jets, Cross-Section Ratios, and the Strong Coupling 

Absolute jet cross sections as used in the above CMS example are subject to many experimental sources 
of uncertainty, and in particular to details of the jet energy calibration. Since the jet-px spectrum is 
steeply falling (see figure [ 6 ]), an uncertainty in the jet energy scale (JES) of 5% would translate into an 
uncertainty on the measured cross section of about 25-30%. A means to at least partially cancel such 
experimental uncertainties was investigated by ATLAS in the form of the ratio of inclusive jet cross 
sections measured for two different center-of-mass energies (2.76 and 7 TeV) [2S]. Figure [ 8 ] shows the 
considerably reduced experimental uncertainties of this ratio compared to the cross sections themselves, 
thus demonstrating the possibility to provide stronger constraints on PDFs. For comparison, the PDF 
uncertainties for various PDF sets are indicated as well. 

Finally, the possibility to extract the strong coupling constant a s (Mz) from inclusive jet production 
data is studied. As can be seen from figure [7] and the DGLAP equations the appearance of gluons 
and the strength of a s (Q ) are coupled. Therefore, a simultaneous fit of PDFs and a s (Mz) with DIS data 
alone is not possible, which is also the reason why the value of a s (Mz ) is fixed in the above-mentioned 
PDF fit of Ref. [1051 . The inclusion of jet data facilitates such a simultaneous fit of PDFs and the 
strong coupling as demonstrated by CMS EDI- In the following, the focus shall be, however, on the 
running behaviour of a s (Q) according to equation ([3]). To eliminate or reduce numerous experimental 
sources of uncertainty—such as those related to the luminosity measurement or the jet energy scale— 
again a ratio is used: that of the inclusive 3-jet over the inclusive 2-jet production cross section, R 32 . 
In this quantity, also theoretical uncertainties like those stemming from the PDFs are reduced. Since 
3-jet production requires one real emission more than the usual dijets, the numerator is approximately 
proportional to a^, while the denominator is oc a 2 S . R 32 is therefore a direct measure for the strength 
of the strong coupling. The CMS collaboration has investigated this ratio with the result a s (Mz) = 
0.1148 ± 0.0014 (exp) ± 0.0018(PDF) ± 0.0050(theo) ill9| . By performing the analysis separately for 
different regions of the average px of the two leading jets, the running could be tested up to scales 
of 1.39 TeV. Figure |9](a) shows the ratio i ? 32 measured by CMS in comparison to predictions at NLO 
obtained with the CT10 PDF set for a series of assumptions on a s (Mz). Figure[9](b) presents a summary 
of extractions of a s (Q) including the latest results achieved at the LHC with scales Q reaching beyond 
1 TeV. 
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Fig. 9. (a) R 32 as measured by CMS in comparison to predictions at NLO with the CT10 PDF for a series 
of assumptions on a s (Mz). (b) Summary of extractions of a s (Q) from low to very high scales Q from various 
experiments including the latest results achieved at the LHC up to Q « 1 TeV. (Adapted from Refs. 11001 1191 
and auxiliary material.) 


6 NNLO: The Quest for Precision 

For a few benchmark processes and observables, the experimental accuracy of QCD measurements at 
the LHC is aiming for the few-percent level or even better. Examples are inclusive single-jet or dijet 
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production, vector boson (pair) production, and top-quark or Higgs-boson production. Clearly, for these 
channels theoretical predictions at NNLO in the strong coupling are needed, as otherwise the extraction 
of Standard Model parameters will be limited by theoretical uncertainties. 

The evaluation of NNLO corrections for a given process requires the calculation of three contributions: 
i) generic two-loop graphs; ii) one-loop graphs with one extra final-state parton; iii) real-emission tree 
diagrams with two additional partons. All of these contributions exhibit infrared singularities, and only 
their sum yields finite results. In explicit calculations, the divergences are typically treated in dimen¬ 
sional regularisation. Similar to what was discussed for the treatment of infrared singularities appearing 
in NLO corrections before, there exist now several subtraction schemes to handle these divergences con¬ 
sistently (12011123) . 

Concerning the above-mentioned processes in need of precision calculations, full NNLO results exist 
for inclusive W/Z -boson production (124111~26j . for cliboson (127)1129) and diphoton production [1301 . 
for top-quark pair production 131 1 and for various Higgs-boson production channels (13214139) . For the 
NNLO evaluation of dijet production, results have been achieved for two-gluon final states, initiated by 
gluon annihilation m, and for quark-antiquark initial states mu. For the inclusive jet cross section 
at yfs = 8 TeV, considering anti-fcj jets with R = 0.7, transverse jet momenta of pr > 80 GeV and 
jet rapidities |y| < 4.4, the authors find an increase of the NNLO prediction with respect to the NLO 
estimate of about 27-16%. The corrections are largest for small jet transverse momenta and decrease for 
larger values of jet pr. The given calculation relies on the MSTW08-NNLO PDF set [106] and uses a 
dynamical scale choice of /i = fiR = p,F equal to the transverse momentum of the leading jet. The size 
of these corrections highlights the importance of higher-order calculations for LHC precision observables. 
In what follows, two standard “candles” (i.e. well-measured processes or processes with theoretically 
well-understood cross sections) at the LHC—the Drell-Yan processes and diphoton production—shall be 
discussed in some more detail. 


6.1 Inclusive Vector-Boson Production 

Due to the relatively simple colour-singlet production mechanism and an experimental signature based 
on leptons that can be measured very accurately, the production of W and Z bosons is well suited for 
precision tests of QCD and electroweak interactions. 

At hadron colliders, massive electroweak bosons are dominantly produced via quark-antiquark anni¬ 
hilation. The production of Z bosons is commonly referred to as Drell Yan process, since in 1970 S. Drell 
and T.M. Yan applied the developing parton model to predict the scaling of the pp —>■ *y/Z —>• p, + p~ 
cross sections da/dM 2 oc I\ L(M 2 /s) 1 with M the invariant boson mass, L the parton luminosity and 
K a constant that depends on the parton spin and coupling nn. In modern QCD, the cross section can 
be expressed as da/dM 2 = (dna^ M ) / (3n c AI 2 s)L qq (M 2 / s), where L qq is the quark-antiquark luminosity, 
«em the electromagnetic coupling, and n c the number of colours. 

Vector-boson processes are relatively simple to calculate, since there is no QCD interaction involv¬ 
ing the final-state particles. In fact, the Drell Yan process was the first hadron-collider channel known 
to NNLO accuracy |124H126j . The NNLO Drell-Yan production cross sections can be calculated with a 
precision of 1% for the renormalisation-scale and factorisation-scale uncertainty. The calculations also 
include the decay of W and Z bosons to leptons. For inclusive Z-boson production at the LHC, the 
NLO-to-LO correction is sizeable (about 20%), while the NNLO-to-NLO correction amounts to 2% only 
. Two public codes implement the NNLO QCD calculations for inclusive Z-boson and VF-boson produc¬ 
tion: Dynnlo (126) and Fewz (14311144] , Furthermore, the dominant electroweak corrections have been 
determined (14514150) . 

The precise measurement of the inclusive IV-boson and Z-boson production cross sections at LHC 
and their successful comparison to precise QCD calculations was one of the first and very important 
confirmations of pQCD at very high energies and momentum transfers. Measurements at fs = 7 TeV 
have been performed with the 2010 data set corresponding to an integrated luminosity of 36 pb _1 by 
ATLAS m and CMS (152) . CMS also measured the cross sections at yfs = 8 TeV using a data set of 

18.2 ft" 1 pm 

The cross sections are measured by counting the number of events in the detector acceptance and 
subtracting the estimated background contributions. The efficiency e is estimated from Monte Carlo 
simulations and corrected for differences between these simulations and the data. To better separate 
experimental and modelling uncertainties (in particular from the extrapolation of the detector acceptance 
to the total phase space), the efficiency is decomposed in a fiducial acceptance, A, and a correction for 
detector effects, C , i.e. e = A ■ C. 
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W+ -A l + v 

W~ -A l~u 

z -a tr 

7 TeV 

ATLAS 

6.05 ±0.12 ±0.21 nb 

4.16 ±0.10 ±0.14 nb 

0.937 ±0.02 ±0.03 nb 

CMS 

6.04 ±0.10 ±0.24 nb 

4.26 ±0.08 ±0.17 nb 

0.974 ± 0.02 ± 0.04 nb 

NNLO QCD 

5.98 ±0.3 nb 

4.2 ± 0.2 nb 

0.991 ±0.05 nb 

8 TeV 

CMS 

7.11 ±0.14 ±0.18 nb 

5.09 ±0.12 ±0.13 nb 

1.15 ±0.02 ±0.03 nb 

NNLO QCD 

7.12 ±0.2 nb 

5.06 ±0.13 nb 

1.13 ±0.04 nb 


Table 1 . Summary of the cross-section results of the inclusive W-boson and Z-boson analyses of ATLAS and 
CMS. The total uncertainty (statistical, experimental and acceptance modelling added in quadrature) and the 
luminosity uncertainties are also given. The theory cross sections are calculated with Fewz (l43i 11441 with the 
MSTW2008-NNLO PDF set 1106 1 and contain renormalisation-scale and factorisation-scale uncertainties and 
PDF uncertainties. 


The fiducial acceptance A is given by the ratio of the number of events passing the kinematic cuts 
applied on particle-level (i.e. without detector simulation) over the total number of generated events. 
The kinematic selection cuts are set close to the requirements on the reconstructed objects typically for 
leptons: px > 20 GeV and \r)\ < 2.5. The detector correction factor C is defined as the number of selected 
events in the sample including the detector simulation over the number of events passing the fiducial 
acceptance A. Typically A is 0.45-0.50 and C is 0.7-0.8, depending on the analysis. 

A summary of the inclusive Drell-Yan cross-section measurements is shown in table [lj The statistical 
uncertainty amounts to about 0.3% for TT-boson production and 0.7% for Z-boson production. The 
systematic uncertainty is dominated by the knowledge of the integrated luminosity, which at yfs = 
7 TeV results in an uncertainty of 3.4% (ATLAS) and 4.0% (CMS) and in 2.5% at y/s = 8 TeV. The 
measurements of ATLAS and CMS agree within their uncertainties, and both agree with the NNLO 
prediction that has an uncertainty of about 5% (dominated by the PDF uncertainty). 


6.2 Differential Vector-Boson Cross Sections 

Apart from the total W and Z production cross sections, also some differential cross sections are known to 
order 0(a 2 )\ this includes the dilepton invariant-mass spectrum, the vector-boson’s transverse momentum 
and its rapidity distribution. 

As an example, figure [To| presents a CMS measurement of the Z-boson production cross section as a 
function of the invariant mass M of the two leptons from the Z decay at yfs = 7 TeV jl54j . The cross 
section falls over eight orders of magnitudes in the mass range of 15 < M < 1500 GeV and clearly shows 
the Z-boson resonance at 90 GeV. The NNLO calculation using the CT10 PDF set describes the data 
well. Similar measurements were performed by ATLAS H55J. 

The boson’s transverse momentum is an interesting observable to test various regimes in the strong- 
interaction dynamics. The high-px part of the spectrum can be modelled by perturbative QCD calcula¬ 
tions considering processes with additional jets ( W/Z + n-jet processes). Presently calculations can be 
performed with up to n = 5 additional partons. When the transverse momentum of the vector boson is 
much smaller than its invariant mass reconstructed from the two leptons, soft QCD radiation is enhanced 
and fixed-order perturbation theory has to be supplemented by the resummation of large logarithmic 
corrections of the form a" In m (px/M) to all orders of a s . The resummed cross sections also include a 
non-perturbative component at momentum scales below 1 GeV. This region can also be modelled using 
Monte Carlo generators implementing parton showers and models for hadronisation. 

To lowest order, W and Z bosons are produced via quark-antiquark annihilation, i.e. qq -A Z. Indeed, 
this process dominates for inclusive Z production within a rapidity range of \y\ < 2.1. However, already 
for Z transverse momenta around 25 GeV quark-gluon scattering qg -A Zq is of similar size, and around 
Pt = 100 GeV the latter process constitutes 80% of the total cross section m- For Z-boson transverse 
momenta of 180 GeV, parton-momentum fractions Xi/ 2 = (M/y/s) exp(±p) of about x = 0.05 are probed. 

The program ResBos |158| implements soft-gluon resummation at low px at next-to-next-to-leading- 
logarithm (NNLL) accuracy and matches to a matrix-element calculation of order 0(a 2 ) at high px- It 
provides fully differential cross sections as a function of the rapidity, the invariant mass and the transverse 
momentum of the vector boson. It also attempts to estimate non-perturbative corrections I159j . Recently, 
first calculations for Drell-Yan lepton-pair production at NNLO accuracy including parton-shower effects 
have been presented [160[ 1161 ]. 
















Carli, Rabbertz, Schumann: Studies of Quantum Chromodynamics at the LHC 


17 


CMS 




' 1 1 1 I T T T T I ' ' T T 1 

-(b) LHCb, 7 TeV 


i 1 1 t 



1 data s , 

t ° MSTW08 

- 

0.4 


data IO 

v NNPDF23 

— 




* CT10 

- 

0.2 

■A V.’bR. Q.?, 


□ ABM 12 

x 

0 



♦ HERA 15 

: 



a JR09 

x 

-0.2 




a 



"T" 


-0.4 

p T > 20 GeV/c 



x 




f ™i ^ 

i- 

0.05 




- 

-0.05 

' ■ . _ 1 _ . _ 1 ■ ■ _ 1 _ 1 _l_i_ 1 _ 1 _ 1 _ 1 _ 

M 1 

■^ 7 ,} ? 

f 


>1 


Fig. 10. (a) Drell Yan cross sections at i/s = 7 TeV as a function of the dilepton invariant mass M. (b) Lb-boson 
charge asymmetry as a function of the lepton pseudo-rapidity. Data are compared to NNLO QCD predictions 
using various PDF sets. (Adapted from Refs. IJ54\UM-) 




Fig. 11. Normalised transverse momentum of the Z boson measured within the detector acceptance at i/s = 
7 TeV. (a) ATLAS data compared to NNLO QCD calculations based on Dynnlo [126] and Fewz [14311144] ! with 
different renormalisation-scale and factorisation-scale settings, (b) Ratio of an NNLO QCD prediction including 
NNLL resummation to the data. Shown are calculations that differ in the parametrisation of non-perturbative 
effects. The bands indicate the respective theory uncertainty estimate. The band around 1 denotes the data 
measurement uncertainty. (Adapted from Ref. '( 16Sb if J 


Early measurements of differential cross sections at ^fs = 7 TeV for Z [163] and W [164 1 153] bosons 
were based on an integrated luminosity of 35-40 pb -1 . Recently, the ATLAS experiment updated this 
measurement to a larger data set corresponding to an integrated luminosity of 4.7 fb -1 . This measure¬ 
ment reaches Z-boson transverse momenta of 800 GeV for rapidities up to 2.4 m- The measurement 
uncertainty amounts to 1% for px < 150 GeV and rises to about 5% for the highest px- The leptonic 
decays of the Z boson to electrons or muons are used in the analysis, and electrons or muons with 
px > 20 GeV and rapidity |y| < 2.4 are required. The invariant mass of the clilepton pair must lie be¬ 
tween 66 < M < 116 GeV. Figure flT] shows the comparison of the data to the NNLO QCD predictions 
based on Dynnlo il26] and Fewz (llTfl [1441 using the CT10 PDF set. Shown as a band are theory 
uncertainties that are due to the renormalisation-scale and factorisation-scale variations and associated 
with the choice of PDFs for both calculations. The measurement uncertainties are below a percent in most 
bins and are smaller than the theory uncertainties that reach 10% around 50 GeV. The two calculations 
are in excellent agreement with each other for px < 100 GeV, but differ at large px by about 10%. In the 
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region 10 < px < 100 GeV the predictions are 10% lower than the data. For px < 10 GeV the fixed-order 
calculations diverge and clearly disagree with the data. 

The resummed calculation based on ResBos using the CTEQ6M PDF [l66j describes the data within 
10%. The observed deviations are covered by the theory uncertainties of the standard non-perturbative 
parametrisation (denoted “BLNY”). A recent new parameterisation denoted “GNW” |159| predicts a 
somewhat different shape. 

An NLO QCD calculation based on Powheg interfaced with Pythia for parton showers and hadro- 
nisation describes the data within 5% for px < 60 GeV, but deviates by up to 20% over the full px range. 
An NLO QCD calculation based on MC@NLO interfaced with Herwig for parton shower and hadroni- 
sation deviates by up to 40% at high p T . Leading order Monte Carlo generators based on multi-leg matrix 
elements, like Sherpa or Alpgen, agree with data within 5% except for the highest px bin. A dedicated 
tuning of parton-shower parameters in Pythia8 interfaced to Powheg jl62j achieves agreement within 
2% for px < 50 GeV. Since the dominant VF ± production mechanisms at the LHC are du —> W~ and 
ud — > W + , measurements of LF ± -boson production cross sections provide a powerful tool to study the 
parton density functions of d and u quarks in the proton. In particular, the measurement of the lepton 
charge asymmetry Ae = (aw+^e+v ~ cr w~^i-p)/{ cr w+^i+ v + <J w-^e.-v) gives valuable information on 
the d-to-u ratio and also on the sea-quark and sea-antiquark distributions, in particular of the strange 
quark. By measuring differential cross sections, for instance as a function of the lepton rapidity, different 
parton-momentum fraction values x can be probed, since the parton-momentum fraction depends on the 
vector-boson mass and on the rapidity y: X \/2 = {AIw,z / y/s) exp (±p). 

ATLAS published such differential cross-section measurements da-^± /z / dye for electrons and muons 
up to lepton rapidities of ye = 4.9 based on the 2010 data set at y/s = 7 TeV 11511 . Leptons with a 
px > 20-25 GeV are selected. For Z bosons, the accuracy reaches about 2% in the central region of the 
Z-boson rapidity and 10% at yz = 3.2. A precision of about 2% is obtained for LF-boson cross sections 
measured within \ye\ < 2.5. For the Ae measurement the accuracy ranges between 4 and 8%. The data 
can be described by an NNLO calculation using the HERAPDF1.5 and the MSTW08-NNLO PDF sets. 
Further measurements by ATLAS and CMS can be found in Refs. 11521 165 . il'67[{1701 . 

The LHCb experiment measured TV-boson production in the forward region covering lepton rapidities 
in the range 2 < \ye\ < 4.5. A measurement was performed with the 2011 data set at yfs = 7 TeV using 
an integrated luminosity of 1 fb _1 . Figure [To^ b) shows the measured lepton charge asymmetry for muons 
with px > 20 GeV as a function of the muon rapidity compared to NNLO calculations using various 
PDF sets. The band shows the statistical (inner, dark) and total uncertainty (light) on the data. The 
experimental accuracy is similar to the theoretical uncertainty estimate, and the predictions are in good 
agreement with the measurement. 


6.3 Production of Photon Pairs with Large Invariant Mass 

At leading order in pQCD, the production of photon pairs can rather easily be described by quark- 
antiquark annihilation, (gg —» 77 ), as for the Drell-Yan process before. However, in contrast to the 
massive vector bosons W and Z , which are predominantly produced in hard scattering processes and can 
easily be identified via their leptonic decay modes, massless photons can also be radiated off initial-state 
or final-state partons or may be created in decays. In particular the 7 r° and p mesons, and to some 
extent also the w, all of which are copiously produced within jets, have neutral decay modes. Since these 
particles are boosted within the jets, their decay photons are collimated and difficult to distinguish from 
a single highly energetic photon. Without applying further selection criteria, the mode of “photon” pair 
production at large invariant mass via meson decays within two (different) jets is even dominating by 
several orders of magnitude. 

To distinguish the “non-prompt” decay photons within jets from “prompt” photons stemming from 
the hard scattering, isolation criteria are applied. Typically, isolation demands that the transverse energy 
not associated to a photon candidate but deposited within a cone of AR = \J ( Ar\) 2 + (Aej)) 2 = 0.4 around 
it may not exceed a few GeV. In addition, a minimal separation of the two photons is required to avoid 
overlap: Z\i ? 77 > 0.40, 0.45 for ATLAS and CMS, respectively. In practice, the experimental selections 
are more complicated (see e.g. Refs. [l7llH74j L and great care has to be taken to ensure that they can 
be matched properly to selection criteria applicable in theory (see Ref. m for a recent discussion on 
photon isolation). 

Finely segmented electromagnetic calorimeters are of great advantage in such photon studies. Inves¬ 
tigating the shower shape of energy deposits in these calorimeters helps to differentiate between elec¬ 
tromagnetic showers that are caused by a single high-energetic signal photon and those from multiple 
photons like in n° —> 77 decays. Photon candidates are separated from electron (or positron) candidates 
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by using reconstructed tracks. Depending on the amount of material to be traversed (e.g. the beampipe 
or the silicon-pixel and silicon-strip detectors), however, highly energetic photons may convert into e^e - 
pairs before reaching the electromagnetic calorimeter. Dedicated reconstruction methods are applied to 
avoid losing these converted photons, which can make up to 50% [ 174] of the total photon yield. 
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Fig. 12. Photon-pair production cross section at 7 TeV centre-of-mass energy as a function of (a) the diphoton 
transverse momentum px, 77 and (b) the azimuthal angular separation Z\^> 77 as measured by ATLAS |172: . The 
data are compared to parton-shower-improved LO predictions by Pythia6 and Sherpa. The MC predictions are 
rescaled by a factor of 1.2 to match the total cross section observed in data. (Adapted from Ref. jn%.) 


For background subtraction, techniques relying on data are preferred. Detailed studies, necessary to 
further quantify detector effects on the photon reconstruction and isolation and to unfold the signal yields, 
are performed with simulated events provided by the LO event generator Pythia {51], 176] or by multi¬ 
jet-improved event generators that include additional real emissions (extra jets) like Sherpa [76) Il77j or 
MadGraph m- The latter is combined with Pythia for parton showering, hadronisation, and for the 
modelling of the underlying event, see section [9} 

At next-to-leading order, the gluon-initiated processes ( gq(q ) —► 77 q{q)) join the annihilation process 
(i qq —> 773 ) together with corresponding virtual corrections. Here, a further complication arises from the 
collinear fragmentation of a hard outgoing parton into a photon. Theoretically, this process is described 
using fragmentation functions D q ^(z, /j,f) and D g _ > 7 (z, gp)- In the limit of gp ~ Q , where Q is the 
scale of the hard scattering, this fragmentation process results in a contribution to the cross section that 
is similar in size to the one of the LO annihilation process D3SI- 

Prompt photons not emerging via fragmentation from the hard process are also called “direct”. The 
effect of the above-mentioned isolation conditions, which aim to suppress fragmentation photons in favour 
of direct ones, have to be properly estimated not only experimentally for mesons decaying within jets, 
but also theoretically, e.g. via FFs in the perturbative calculations. The NLO parton-level program Res- 
Bos {17911180] effectively includes the fragmentation of one quark/gluon to a single photon at leading 
order and additionally features soft and collinear gluon resummation (see section [8]). Diphox [T78I pro¬ 
vides parton-level results at NLO taking fully into account up to two fragmentation photons. Although 
formally an NNLO box graph, the process (gg —► 77 ) is drastically enhanced at the LHC through the 
large gluon luminosity. It is comparable in size to the LO terms, and therefore Diphox includes this con¬ 
tribution at NLO precision, i.e. up to N 3 LO corrections in the strong coupling a Sl via Gamma2mc m- 
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Finally, a full NNLO calculation is available in the form of the 2yNNLO [130] program, however without 
consideration of fragmentation photons. 


CMS <s = 7 TeV 5.0 fb' 1 


CMS Vs = 7 TeV 5.0 fb' ! 



Atp.,, [radian] 


A<p yy [radian] 


Fig. 13. Photon-pair production cross section at 7 TeV centre-of-mass energy as a function of the diphoton 
azimuthal angular separation Af> 71 as measured by CMS EH- (a) Comparison of the data to Sherpa using 
tree-level matrix elements matched to the parton shower and to NLO results from ResBos that also include 
soft-gluon resummation, (b) Comparison to full NLO calculations including contributions of photons produced 
in jet fragmentation from Diphox + Gamma2mc as well as to NNLO results from 2yNNLO. (Adapted from 
Ref. fTffl .) 


Four kinematic variables are usually chosen to investigate the differential photon-pair production 
cross section: the invariant mass m 77 , the transverse momentum of the photon pair Pt^j, the azimuthal 
angular separation Z\</> 77 , and the cosine of the polar angle in the Collins-Soper reference frame [182] of the 
diphoton system, cos0* 7 . The invariant mass is of obvious interest for resonance searches, where cos0* 7 
is useful to examine the spin of diphoton resonances. The transverse momentum p t i77 and the quantity 
A(f) 77 are well-suited for the comparison of specific aspects of the theoretically very challenging description 
of photon-pair production. In particular, at low PT,~a or for well-balanced photons at A<f> 77 ss 7r, where 
multiple soft-gluon emission becomes important, fixed-order calculations are not expected to work. As 
an example, figure [T2] shows a comparison between data and parton-shower improved LO predictions 
by Pythia6 and Sherpa for PT,-yy and A<j> 77 from ATLAS [172] . Since both MC event generators fail 
to accurately predict the total cross section, the distributions have been rescaled by a factor of 1.2 to 












































Carli, Rabbertz, Schumann: Studies of Quantum Chromodynamics at the LHC 


21 


match the total cross section in data. In both cases Sherpa provides a good description of the whole 
distribution, while Pythia6 exhibits some discrepancy at low PT,-yy 


Figure 13 ’a) furthermore presents comparisons for the cross section differential in Z\<(> 77 for Sherpa 


and to NFC) from ResBos including soft-gluon resummation. Again Sherpa gives a good overall de¬ 
scription, while ResBos predictions are too low by factors of 2-3 at small separation angles. A similar 
behaviour can be seen in the comparison to NLO from Diphox + Gamma2mc in figure [l3[b), although 
due to resummation ResBos performs somewhat better around Z\0 77 ss 7r. The best description of the 
data is given by the NNLO result from 2yNNLO (figure |l3|(b)), even though the fragmentation contribu¬ 
tion is not included. ATLAS observes similar results when comparing their data to the predictions from 
Diphox + Gamma2mc and 2yNNLO [1721 . 
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7 Multi-Jets: Precision Meets Multiplicity 

Understanding the QCD dynamics governing the emergence of multi-jet final states is of utmost impor¬ 
tance at the LHC. Given that typically jet objects with minimal values of their transverse momentum 
of order 30 GeV are considered, the phase space for producing multiple jets is huge. Accordingly, these 
jet-production processes need to be described as well as possible, as they constitute severe backgrounds 
to almost every search for new physics and furthermore have an impact on the appearance of the signals 
themselves. The modelling of inclusive processes that receive contributions from final states with poten¬ 
tially many jets is the realm of so-called “matching and merging” techniques that combine (multi-parton) 
fixed-order calculations with parton-slrower simulations. 

On the one hand, fixed-order calculations provide a well-defined estimate for inclusive production cross 
sections and reliably account for hard, well-separated parton configurations. On the other hand, parton- 
shower simulations capture the dominant terms driving the emission of additional soft and/or collinear 
partons and thus provide the means to perturbatively account for the internal structure of jets and 
inter-jet energy flows. When combining fixed-order matrix elements with parton-shower resummation 
beyond leading-order 2-to-2 scattering processes, the obstacle of double-counting configurations that 
appear in both approaches needs to be resolved. Furthermore, any consistent scheme should neither spoil 
the logarithmic accuracy of the inherent parton-shower resummation nor destroy the formal precision of 
the fixed-order part. 

There are two basic strategies to distinguish. Tree-level merging techniques correct the hardest emis¬ 
sions of the parton shower off a given core process through exact leading-order QCD matrix elements 11831 - 
1185) . This is achieved through a slicing of the real-emission phase space in terms of a hardness-measure 
that regulates any infrared singularities and allows to consistently combine matrix elements of vary¬ 
ing parton multiplicity dressed with parton showers into an inclusive sample. Several variants of such 
leading-order merging techniques exist and are widely used in LHC analyses, see Refs. [6li Il86lfl89j . 

The second ansatz relies on the exact matching of an NLO QCD calculation with a parton-shower 
cascade off the underlying Born process. Accordingly, the real-emission correction as part of the fixed-order 
calculation has to be properly synchronised with the first, i.e. hardest, shower splitting. Furthermore, the 
NLO accuracy with respect to the inclusive production process considered needs to be preserved. Two 
basic solutions exist to this problem, which are known as MC@NLO EM] and Powheg m- Over the 
last years there has been tremendous development in implementing these techniques for a wide range of 
processes and ultimately their automation, see Refs. [TT] 11921 - 1195) . 

Most recently hybrid solutions emerged that combine NLO plus parton-shower calculations with 
higher-order tree-level QCD corrections 1196 , 197) or that even combine parton-shower-matched NLO 
calculations of varying jet multiplicity 119511202] . 

In the following section, the main focus shall be on a class of processes that constitutes a prototypical 
Standard Model background and serves as a test bed for the various types of advanced QCD calculations 
outlined above: W/Z+ jets production. 


7.1 Weak Bosons and Jets 

In the production of W and Z bosons in association with jets (LF+jets), very high jet multiplicities can 
be reached at the LHC. Cross sections for the associated production with up to seven jets with transverse 
momentum px larger than 30 GeV have been measured by ATLAS and CMS. Already the leading-order 
prediction for this process is of order a 7 s . However, this class of processes is extremely important as it 
constitutes a major background to many other processes with complex final states such as top-quark pair 
or diboson production or rare signatures from the Standard Model or beyond. 

The LF+jets processes allow an assessment of the validity of Monte Carlo generators like Alpgen [6j, 
Sherpa jTSJITB], or MadGraph pro] (together with recent progress in MadGraph5_aMC@NLO ITT] ) 
that merge tree-level matrix-element calculations with up to five additional partons interfaced with parton 
showers and hadronisation. Presently, many experimental analyses are in a transition phase adopting NLO 
matrix-element calculations plus parton showers as baseline for MC event generation instead of the former 
LO based ones (see Section |5j. 

Figure [l4| shows the VF+jets cross section as a function of the leading-jet px as measured by CMS [203) 
and ATLAs |204| using proton-proton collision data recorded at y/s = 7 TeV. In the ATLAS measure¬ 
ment, jets are defined by the anti-fc t jet algorithm using a radius parameter of R = 0.4 and measured 
within 30 < px < 300 GeV and \y\ < 4.4. The results for up to four additional jets are superimposed and 
are compared to the predictions of Sherpa and Alpgen. The measured shape of the jet px distribution 
is described by the Monte Carlo simulation within the uncertainties. 
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Fig. 14. The W +jet cross section as a function of the transverse momentum of the leading jet as measured by 
CMS (a) and ATLAS (b). The data are compared to various QCD calculations. For the ATLAS measurement, jet 
multiplicities up to four jets are shown. The results for higher jet multiplicities are scaled by the factors given in 
the figure. The leading-order multi-leg calculations are normalised to the inclusive NNLO cross section. (Adapted 
from Refs. 12031 \204l -) 


In the CMS measurement, jets defined by the anti-fc* jet algorithm with R = 0.5 are selected with 
Pt > 30 GeV and |y| < 2.4. The jet spectrum reaches up to pt values of 800 GeV. The data fall below 
the leading-order QCD predictions towards high p x- The MadGraph (Sherpa) estimate using up to 
four additional partons is almost 50% (20%) higher than the data for px > 200 GeV. 

At the time of the LHC start-up in 2008, NLO QCD calculations with maximally two associated jets 
were available 2T)51 . Since then, in particular, with the BlackHat-Sherpa [90) ' 2061 program significant 
progress was made in NLO QCD calculations for W or Z production in association with up to five 
additional partons ISB- Fixed-order calculations, by nature, are not matched to parton showers and 
hadronisation, so that corrections to account for these effects are derived from generators like Pythia 
or Herwig are needed. The resulting predictions from BlackHat-Sherpa are in good agreement with 
the measured distributions within the systematic uncertainties as shown in figure |14| 

Recently, automated computations of NLO QCD cross sections matched to parton-shower simulations 
for a large variety of processes (with up to four partons) have been developed. The programs Sherpa 
(in the Meps@NLO option [8, l98j ) and MadGraph5^\.MC@NLO [74] presently provide the most 
accurate predictions for multi-jet final states at particle level. 

Similar measurements 21)7! for the Z+jets processes using the full data set at y/s = 7 TeV are 
presented in figure [T5| as a function of the leading-jet px and the fourth-leading jet px- Jets are defined 
using the anti-fc t jet algorithm with R = 0.4 and are selected if px > 30 GeV and \y\ < 4.4. The cross 
sections are normalised to the inclusive Z cross section. Again, the BlackHat-Sherpa predictions are 
in good agreement with both measured distributions. 

An NLO QCD calculation for the inclusive Z cross section merged with parton showers and hadroni¬ 
sation as incorporated in MC@NLO is capable of describing the data for events with one or two low-px 
jets, but falls far below the data towards high values of the leading-jet px- This indicates that the fraction 
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Fig. 15. Z+jets production cross section as a function of the leading and fourth-leading jet transverse momentum 
as measured by ATLAS. The data are compared to various QCD calculations. (Adapted from Ref. \207j.) 


of events exhibiting a second jet increases considerably with the px of the leading jet and that the parton 
shower approximation, which is used by MC@NLO for the second jet, fails to describe the data. 

The results are also compared to two multi-leg Monte Carlo simulations produced with up to five 
additional partons. The Alpgen simulation overestimates the data at high px for the leading jets and 
underestimates the data for the fourth-leading jets. Sherpa is in good agreement with the data con¬ 
sidering the experimental uncertainties. The measurements of the Z+jets cross section are limited by 
the experimental uncertainty on the jet energy measurements. By considering jet-multiplicity ratios, 
many experimental and theoretical uncertainties cancel, allowing for an even more precise data-to-theory 
comparison. 

The ratio of the Z+ jet cross sections for iVj et + 1 over 7Vj et (R(n+i)/n = cr(Z + lVj et + 1 )/a(Z + iVj et )) 
was measured by ATLAS [W] , An earlier measurement was performed by the CMS experiment on 
a smaller data set [2081 Figure [l6| shows the ratio R(n+i)/n for successive exclusive jet multiplicities 
for all events (“inclusive”) and for events with at least one jet with px > 150 GeV. For the inclusive 
case, the ratio shows a rather flat, approximately linear dependence on the jet multiplicity with a small 
slope only, while for events with px > 150 GeV the measured distribution steeply rises towards low 
jet multiplicities. These measurements illustrate two limiting cases of scaling patterns, the “staircase 
scaling” |209j with R(n+i)/n constant and the “Poisson scaling” [210J with R(n+i)/n = (N)/N, where 
( N ) is the average number of jets. These scaling patterns can be qualitatively understood from the 
expected Poisson-distributed jet multiplicity. For large jet-px differences, the Poisson-scaling is directly 
observed (see figure [l6](b)). However, for low px and small _/Vj et , the emission of additional partons is 
suppressed by the PDFs, and for high A( et the emission of additional partons no longer follows a Poisson 
distribution due to the non-Abelian nature of QCD, leading to a proliferation of the number of jets 
originating from gluon splitting (see figure [l6|( a)). 

The scaling patterns observed in data are described by the Monte Carlo simulations BlackHat- 
Sherpa, Alpgen and Sherpa. The scale uncertainty is only shown for the NLO QCD calculation based 
on BlackHat-Sherpa. The different grey shades indicate the scale uncertainty when it is considered 
uncorrelated (as proposed in Ref. [211j f or correlated across the various jet multiplicities. Assuming the 
predicted and now observed scaling patterns to be valid, it becomes possible to give estimates for even 
higher jet multiplicities where there exists no complete NLO prediction yet [212] . 
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Fig. 16. Ratios of Z+jets cross sections for successive exclusive jet multiplicities for (a) inclusive events and (b) 
for events with at least one jet with px > 150 GeV and |j/j e t| < 4.4. The data are compared to various QCD 
calculations. (Adapted from Ref. \207(.) 


7.2 Weak Bosons and Jets with Flavour 


Several key processes in the Standard Model and beyond—such as the production of single top quarks, top- 
quark pairs or Higgs bosons—involve heavy quarks. The understanding of their production mechanisms 
is therefore an important part of the LHC physics programme. For a detailed discussion of flavour physics 
at the LHC see Chap. 8 of Ref. [T] . 

The predictions of processes with heavy-flavour content are more difficult to accomplish than those for 
the inclusive case. Bottom (5) quarks are significantly heavier than the proton, and in strong interactions 
at high scales they can only be created in pairs. The calculations can be classified in two types: In the 
4-flavour (4F) scheme heavy quarks appear only in the final state and are typically considered massive, 
while 5-flavour (5F) scheme calculations include heavy quarks in the initial state, i.e. as present in the 
proton. For this purpose, c-quark and 5-quark PDFs are used. 

The production of W and Z bosons in association with heavy flavour is an important background for 
many processes. Several NLO QCD calculations are available |214d224] . Recently, the W + 55 and Z + bb 
processes (with and without additional jets) have been included in automated NLO calculations matched 
with parton showers and hadronisation m- 


Z bosons provide a particularly clean experimental signature that allows for precise measurements. 
Figure [IT] presents a recent ATLAS measurement of the cross section of Z bosons produced in association 
with 5 jets as a function of the px of the leading jet containing a B hadron. The measurements uses the 
full 2011 data set. Transverse jet momenta up to 500 GeV are reached. The data are compared to various 
QCD calculations. 

The Mcfm prediction is based on a 5F NLO QCD calculation of Z + 15 jet 1217 . 2191 . corrected 
for hadronisation and effects of multi-parton interactions (see section 9.4). Full particle-level predictions 
with NLO QCD matrix element calculations are also obtained using MadGraph5_aMC@NLO [74j in 
both the 4F and 5F schemes. For these calculations, the renormalisation and factorisation scales are set 
to the transverse Z-boson mass and varied up and down independently by a factor of two to assess the 
residual dependence on this scale choice. The data are described within the experimental and theoretical 
uncertainties. 
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Fig. 17. Production cross sections for Z bosons produced in association with jets with heavy-flavour content as 
a function of the jet pr . Superimposed are various QCD calculations. (Adapted from Ref. f213f .) 


Also shown are the predictions for multi-leg LO QCD calculations based on Alpgen using the 4F 
scheme with up to five partons in the matrix element and Sherpa using the 5F scheme with up to 
four partons. The shape of these predictions is in good agreement with the data. However, the LO 
normalisation of these predictions undershoots the measurement, Sherpa reproducing the normalisation 
better than Alpgen. 

Further measurements by ATLAS and CMS of vector-boson production in association with heavy- 
flavour jets can be found in Refs. [225H227] . 


8 Resummation: The Realm of Large Logarithms 

When discussing the process of vector-boson production before, the necessity to incorporate effects of 
multiple emissions of initial-state and final-state partons to appropriately model certain observables has 
been touched upon already. In fact, phase-space regions sensitive to the emission of soft and/or collinear 
QCD quanta often resist a satisfactory description in terms of fixed-order perturbation theory. Rather, 
one is forced to identify the relevant enhanced contributions, namely large logarithmic terms, and to 
reorganise the perturbative expansion for the observable at hand. The aim is to account for the enhanced 
terms to all orders in perturbation theory. This procedure is referred to as “resummation” and allows an 
appropriate treatment in a wide range of perturbative QCD phenomena. 

Examples for hadron-hadron collider observables sensitive to multiple-emission effects include the 
transverse-momentum distribution of lepton pairs or gauge bosons at low p T [1581 12281 1229) . event- 
shape variables [230H232) . or certain jet-substructure observables [233H235] . Furthermore, resummation 
techniques become important when the real-emission phase space is constrained, leaving behind finite 
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but possibly large uncancelled QCD corrections. Examples thereof are observables that exhibit thresh¬ 
olds [236H239] or that rely on an explicit veto on QCD activity [24014244] . Resummed expressions for this 
type of observables might enable us to extract QCD parameters such as the strong coupling, quark masses 
and the parton distribution functions much more accurately than when having to rely on fixed-order per¬ 
turbation theory alone. In the following, two concrete examples shall be considered: the application of a 
veto on additional jet activity and the jet-mass distribution in Z+jet production. 


8.1 Jet Vetos and Gap Fractions 

At first, so called “gap fractions” will be discussed. For this class of observables, used also to select 
the weak-boson-fusion production channel of Higgs bosons (see Chap. 6 of Ref. one considers a 
specific kinematic structure of two hard and widely separated jets. Analysing the QCD activity in- 
between these hard jets, events are vetoed in case there are additional jets with transverse momentum 
above a certain veto scale Q o, filtering out those events that exhibit a gap between the two boundary 
jets. On the one hand, the presence of the potentially soft jet-veto scale Q o induces large logarithms 
with argument Q/Qo , where Q denotes the typical hard process scale, e.g. the boundary jets’ transverse 
momenta. This demands for the all-orders treatment of wide-angle soft-gluon radiation [2401 12411 . On the 
other hand, when considering large rapidity separations Ay between the boundary jets, manifestations 
of BFKL-like dynamics [245[ [246] are expected, requiring to sum terms proportional to a™(Ay) n to 
all orders n 124711248] . In Refs. [249 , 1250 ATLAS presented measurements of dijet production with a 
veto on additional jets. In order to gain further insight into the QCD dynamics determining the gap- 
fraction and gap-jet measurements, Ref. |250j considers additional azimuthal-decorrelation observables. 
In this latter analysis, the data are compared to predictions from Powheg T93; [ 2511 . interfaced to the 
Pythia8 and Herwig event generators, as well as to the program Hej 12'52'j with and without invoking 
the Ariadne shower model [253] . The two Powheg simulations provide NLO accuracy for the inclusive 
dijet production process invoking leading-logarithmic DGLAP resummation through the parton showers 
attached. The Hej approach provides a resummation of small-x, BFKL-type logarithmic terms that can 
be supplemented with DGLAP resummation through, in this case, the Ariadne parton cascade. 

In figure [18] the gap fraction as a function of the veto scale Q 0 is presented for various slices of 
the boundary jets’ rapidity separation Ay. The fraction of events exhibiting a rapidity gap decreases 
as the veto scale is lowered or when increasing Ay. This corresponds to the expectation of an increase 
of the jet activity when considering lower jet transverse-momentum thresholds together with the rise in 
jet multiplicity when considering larger rapidity intervals. The Powheg +Pythia8 prediction yields a 
reasonable description of the data. Employing the Herwig parton shower the agreement worsens, an 
effect that originates from a prediction of too many jets above the veto scale. For the two predictions 
from Hej the agreement with data improves when considering large rapidity separations. Given the high 
precision of the measured data, no single theoretical prediction is capable of describing the data in all 
phase-space regions. Considering more observables, including angular variables, Ref. [250] establishes 
the DGLAP-based Powheg +Pythia8 and BFKL-like Hej +Ariadne simulations to yield the best 
overall description of the data. So far no clear-cut statements about the evidence for BFKL effects 
can be established. However, the small experimental uncertainties allow theoretical models for QCD 
radiation between widely separated and high-transverse momentum jets to be further constrained and 
improved [5TUI . 

8.2 The Jet-Mass Distribution 

Many phenomenological studies over the past years have investigated and highlighted the potential of 
jet-substructure techniques to identify the hadronic decays of boosted heavy objects or to discriminate 
scenarios of new physics from Standard Model (i.e. QCD) backgrounds |258H260j . The successful appli¬ 
cation of jet-substructure methods requires a more detailed and even finer understanding and modelling 
of the inner structure of QCD jets. The new experimental challenges to address include the consideration 
of various jet algorithms and radii that need to be calibrated, the quantification and possible removal of 
non-perturbative contributions from underlying event or pile-up activity, or the application of tracking 
and flavour-tagging methods in busy environments. On the theoretical side, precise and reliable predic¬ 
tions for the new observables at hand are required. This, in particular, asks for resummed calculations, 
in the form of Monte Carlo parton-sliower simulations or as dedicated analytical results. 

ATLAS and CMS have presented several analyses that address these questions and compare their 
data to theoretical models [5544351 12624265] . One particular and rather simple observable shall be con¬ 
sidered as an example here: the invariant mass of jets—the jet-mass distribution. The jet-mass variable 
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Fig. 18. Gap fraction for dijet events as a function of the jet-veto scale Q o for various Ay ranges. The data 
are compared to predictions from Powheg +Pythia8, Powheg +Herwig and from the generator Hej with 
and without inclusion of the Ariadne shower algorithm, (a) Direct comparison with the unfolded data (for 
illustrative purposes the various results are shifted by arbitrary constants); (b) ratio of theory over data. (Adapted 
from auxiliary material provided with Ref. !250/ .) 




Fig. 19. Unfolded jet-mass distribution in Z(££)+ jet events. Jets are reconstructed using the anti-fc t jet algorithm 
with distance parameter R = 0.7 and are groomed using the pruning technique EBB- The data from CMS are 
compared to Monte Carlo predictions from MadGraph +Pythia6 (solid) and Herwig-I — f (dotted). The right 
panel shows the ratio of Monte Carlo over data. (Adapted from Ref. 1262/ . ) 
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probes the composition of jet objects and thus depends on the production process, the jet algorithm and 
its parameters used to cluster the jets, and on additional selection or filtering steps applied to the final 
states. The aim of “jet grooming” procedures is to reduce the impact of underlying-event and pile-up con¬ 
tributions on the jet properties and sometimes also to remove certain types of intra-jet parton splittings, 
all that in order to maximise the sensitivity to the desired signal. Examples of jet-grooming procedures 
include trimming [ 266] . pruning [251] . mass-drop I21TT . I268 J or soft-drop 12651 filtering. Dedicated the¬ 
oretical predictions for the jet-mass distribution exist for a variety of production processes, including 
dijet final states, vector-boson plus jet production and Higgs-boson plus jet production (2701 - 1275) . The 
mass distribution of QCD jets after the application of grooming techniques has recently been discussed 

in Refs. |573112751 . 

In Ref. j2B2 the CMS collaboration has presented an analysis of the jet-mass observable in dijet 
and W/Z+]et events, focusing on the effects of various grooming techniques and comparing data to 
parton-shower Monte Carlo predictions. For an integrated luminosity of 5 fb -1 of 7 TeV collision data 
they considered both the anti -k t and Cambridge-Aachen (276H278) jet algorithms using various distance 
parameters, namely R = 0.7 for anti-fc t and R = 0.8 and R = 1.2 for Cambridge-Aachen. The groomed 
and ungroomed data are compared to Monte Carlo predictions from Pythia6 and Herwig+-|-. As an 
example, the unfolded jet-mass distribution of Z (£(.)+jet events is shown in figure [l9| Depicted there is 
the jet-mass variable for anti -k t jets after a pruning procedure, analysed for four slices in jet transverse 
momentum (note the arbitrary scaling factors). For all the jet-px slices both generators describe the data 
well. A similar level of congruence is observed for the alternative grooming techniques and combinations 
of jet algorithm and distance parameter. Somewhat larger deviations between data and Monte Carlo 
estimates are observed for regions of small jet mass. However, it is in particular this region of phase space 
that is most affected by non-perturbative corrections from pile-up, underlying event and hadronisation. 
Attempts to analytically estimate the effects of the underlying event and parton-to-hadron fragmentation 
on the jet-mass distribution have been presented in Ref. (2791 . 

From the analysis presented by CMS and a similar one presented by ATLAS [561, agreement of jet- 
mass and substructure observables with modern parton-shower Monte Carlo tools is established. This 
paves the way to successfully use these techniques in future analyses and particularly in searches for new 
phenomena. 


9 Beyond Perturbative QCD 

So far the focus has been mainly on measurements that can reliably be addressed using QCD perturbative 
techniques. Experimentally, however, hadrons are observed in the detectors, and accordingly the transition 
of partons into hadrons needs to be accounted for. Furthermore, there is a wealth of effects that require 
the use of non-perturbative methods and models—most prominently the so-called “underlying event” 
that accompanies a hard scattering. As an example, the underlying event includes contributions from 
multiple interactions of the colliding protons’ constituents (“multi-parton interactions”) that give rise to 
the production of many additional final-state particles. 

In this section these aspects shall be exemplified by presenting LHC measurements of jet-shape vari¬ 
ables, of jet-radius ratios and of observables sensitive to soft-particle production and double and multi- 
parton scattering. 


9.1 Jet Shapes 


Focusing on the internal structure of jets, the profile of the distribution of transverse momentum within 
a jet can be examined. For this type of observable the term “jet shape” has been coined |280 i, ; 28l| . The 
differential jet shape p{r) as a function of the distance n = \J\Ai^ et y) 2 + {A it - ]et (t>) 2 to the jet axis is 
defined as the average fraction of jet px contained inside an annulus of inner radius r a = r — Ar /2 and 
outer radius = r + Ar/ 2 for an ensemble of N jets: 


p{r) 


1 1 ^hr a <ri<r b PT,i 

N ^s Ar E n <RPT,i 


where the second sum runs over all jet constituents i. The integrated jet shape < P(r) is then given by the 
integral of the differential jet shape up to a radius r (see the illustrations in figure [20|(b) ). Conventionally, 
measurements are presented in terms of 1 — ^(r = r core ), where r core is usually taken to be 0.3 (282) . 
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Perturbatively, the internal structure of jets is determined by multiple emissions of gluons and depends 
on the type of initiating parton, i.e. a quark or gluon. Perturbative QCD predicts that gluon-initiated jets 
are broader in shape with a higher particle multiplicity on average [2831 . To some extent it is therefore 
possible to differentiate statistically between gluon and quark jets by looking into the jet substructure. 
Jet shapes, however, are also sensitive to non-perturbative effects of the parton fragmentation and the 
underlying event. Corresponding measurements are valuable inputs to improve the modelling of such soft 
contributions in MC event generators like Pythia, Herwig-I —b, or Sherpa. 
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Fig. 20. (a) Differential jet shape p(r) as a function of the distance r from the jet axis for jets with 110 < pt < 
125 GeV from CMS. The data are compared to various tunes of the Pythia6 and to the Pythia8 and Herwig+ + 
event generators, (b) Illustration of the jet shape observables, (c) Integrated jet shape 1 — \P(r = 0.3) as a function 
of the jet pt from ATLAS 12841 . The data are compared to the prediction for gluon-initiated and quark-initiated 
jets using the Pytiiia6 MC event generator with tune Perugia2010. (Adapted from Refs. \28f \ \28S j.) 


Figure |20{a) shows a CMS measurement 12851 of the differential jet shape, for the example of the jet pt 
interval 110 < pt < 125 GeV, in comparison to several LO MC event generators at particle level. While 
Herwig++ and Pythia8 predict somewhat broader jet shapes than measured, the PythiaG tunes 
D6T and Z2 deviate in the opposite direction. The best description of the data is given by Pythia6 
with the tune Perugia2010, something, which is also observed in a similar investigation by the ATLAS 
collaboration |284| . By employing this tune for the integrated jet shape, shown in the form of 1—!A(r = 0.3) 
as a function of the jet p t in figure [20lc), ATLAS demonstrates the sensitivity of this observable with 
respect to the jet-initiating parton. At small jet pt the examined inclusive jet sample is predominantly 
composed of gluon jets, while with increasing jet pt the quark-initiated component grows. In the figure, 
the change in the fraction of quark-initiated to gluon-initiated jets is averaged over the jet rapidity up to 
| y\ < 2.8. If studied double-differentially, a mild dependence on \y\ of this fraction is observed as expected. 
The sensitivity of jet shapes to non-perturbative effects and modelling uncertainties, however, prevents 
an extraction of the quark-gluon jet fraction. Instead, they provide valuable input to the tuning of these 
effects in MC event generators. 

In the context of jets initiated by heavy boosted objects, the investigation of jet substructure has 
attracted much attention lately. Potential new particles like Z' bosons or heavy t' quarks that decay 
lradronically may be at the origin of jets with two or three, respectively, high-energetic partons frag¬ 
menting into only one jet. The required boosts for these hypothetical particles are attainable at the 
LHC, in particular at 13 TeV centre-of-mass energy as foreseen for the restart in 2015. Searches for these 
particles profit from a better understanding of the QCD background in the form of quark-initiated or 
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gluon-initiated jets. For some of the investigated observables, e.g. jet mass, resummed predictions are 
available (see section [8J. 


9.2 Jet-Radius Ratio 


As described in section [5j the inclusive jet cross section can precisely be compared to NLO predictions, 
and the strong coupling constant can be extracted. A small correction for non-perturbative effects, which 
is px-dependent and increases at small px values, has to be taken into account. Looking more closely 
into the results presented by ATLAS for anti-/c t jet sizes of R = 0.4 and 0.6 [3D] or by CMS for anti-fc t 
jet sizes of R = 0.5 and 0.7 [SUES], it is observed that the relative normalisation of the measured cross 
sections to the theoretical predictions exhibits a dependence on R. Theoretically, this R dependence has 
been examined in Refs. [27511255] . where in a collinear approximation it was found that the impact of 
perturbative radiation and of the non-perturbative effects of hadronisation and the underlying event on 
jet transverse momenta scales roughly with Ini?, —1/R, and R 2 for small i?, respectively. By choosing 
the jet size parameter R , one can therefore steer which aspects of jet formation are emphasised in a jet 
analysis. In Ref. [279 i it is suggested, in order to gain insight into the interplay of these effects, to study 
the relative difference between inclusive jet cross sections that emerge from two different jet definitions: 


/ der 
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= 77(alt, ref) — 1. 


Provided that partons in opposite hemispheres are not clustered together (a condition that is usually 
fulfilled), the LO two-parton cross sections are identical for arbitrary jet algorithms. Therefore, only par- 
tonic final states with three or more partons lead to a numerator different from zero. Hence, equation ( |9.2[ ) 
defines a three-jet observable, 77, for which it was shown in Ref. [257j that it is calculable to NLO with 
terms up to of with Nlojet++ :288j 2~89) . 



CMS, 5 fb 1 


VY= 7 TeV 



Fig. 21. Jet-radius ratio 7? as measured by CMS in comparison to (a) fixed-order predictions with and without 
non-perturbative corrections and to (b) particle-level predictions of LO and NLO event generators with matched 
parton showers and modelling of hadronisation and the underlying event at central rapidity |f/| < 0.5. (Adapted 
from auxiliary material provided with Ref. \3$-) 


Studies in that direction have been performed at the HERA collider by the ZEUS collaboration for 
two different jet algorithms [290) and by the ALICE experiment for the two different anti -k± jet sizes of 
R = 0.2 and 0.4 12911 . The results of a CMS analysis for their two jet sizes of 0.5 and 0.7 is presented in 
Ref. El] as a function of jet px and rapidity y. It is expected that QCD radiation reduces this jet-radius 
ratio 7 Z below unity and that this effect disappears with increasing collimation of jets at highpx- Figure [21] 
confirms this assumption with a comparison of the measured ratio 77 as a function of the jet px up to 
the TeV scale to theoretical predictions at central rapidity \y\ < 0.5. Figure [2T[a) clearly demonstrates 
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that fixed-order calculations up to NLO, even when combined with non-perturbative corrections, are 
systematically above the data. The LO event generators Pythia6 and Herwig+-|- lie systematically on 
either side of the data, as can be seen in figure |2l)(b). Presumably, this can be improved by including this 
observable into the tuning of these generators. The best description is given by Powheg that matches 
the dijet production process evaluated at NLO with the parton showers and non-perturbative models of 
Pythia6, emphasising the importance of parton showers for the given choice of jet sizes. 


9.3 Soft Hadron-Hadron Collisions 

When two hadrons collide at high energy, their entire volumes or at least a significant part of them 
participates in the collision. Most of these collisions are caused by so-called “soft” hadronic interactions. 

In a simple geometrical model, the total hadron-hadron collision cross section is related to the trans¬ 
verse extension of the scattering system. Therefore, it is expected to be of the order of cr to t ~ 0(2 fm 2 ) ss 
0(40 mb -1 ). Since the effective range of soft hadron interactions does not change much with energy, only 
a weak energy dependence is expected. Moreover, if the collision is soft, the scattering cross section at 
high energy should not depend on the quark structure of the colliding hadrons. 

In 1952 Heisenberg used a geometrical model where hadrons are modelled as Lorentz-contracted disks 
with a given energy density. The interaction takes place in the disk-overlap region, and only if there 
is sufficient energy to create two mesons. In this model the total cross section increases like In 2 s. A 
consequence of this model is that the interaction cross sections of all hadrons are similar at high energy. 
For a recent review on the model giving a In 2 s high-energy behaviour see Ref. [292] . 

General arguments based on unitarity, analyticity and factorisation imply a bound on the high- 
energy behaviour of total hadronic cross sections. This bound is independent of the details of the strong- 
interaction dynamics and states that the total cross section can not rise faster than In 2 s 12331 . 



Fig. 22. Comparison of total and elastic cross-section measurements in proton-proton and antiproton-proton 
collisions and in cosmic-ray interactions as a function of the centre-of-mass energy. (Adapted from Ref. \29f(.) 


Soft hadron collisions can not be calculated in perturbative QCD. Instead, they are usually described 
by phenomenological models where the force between hadrons is assumed to be mediated by mesons, for 
instance by pions [295]. The basic idea is that the exchanged meson interacts with the colour charges 
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confined in the scattered hadrons. Since many mesons exists, the strong force between hadrons can 
be quite complicated in such models. The typical range of the strong force is anti-proportional to the 
mass of the exchanged meson. Consequently, the largest range is provided by the lightest meson which 
has the appropriate quantum numbers. Such models are usually called “Regge theory”. For reviews see 
Refs. [29611297] . Regge theory aimed at explaining the spectrum of the various hadrons, the forces between 
them, and the high-energy behaviour of the hadron-hadron scattering cross section. It is based on simple 
assumptions like that the scattering matrix is analytic and unitary. Mesons are understood as composite 
quark-antiquark states bound by an effective strong potential that depends on the distance between the 
quark-antiquark pair r and the orbital angular momentum J. Elastic hadron-hadron scattering is viewed 
as the exchange of a meson in the s-channel or f-channel. In 1961, Chew and Frautschi 1298] made the 
observation that for 7r ~p —» n°n scattering there is a linear relation between the mass of the exchanged 
hadron (to) and its angular momentum J. Each integer value of the slope a —called “Regge trajectory”— 
corresponds to one particle that can be experimentally identified (“Regge-pole”). In the time-like region 
where t > 0 (s-channel meson exchange) the momentum transfer t corresponds to the hadron mass to. If 
one continues the linear relation between to and J using the slope a for t < 0, i.e. the space-like regime, 
a good description of the experimental measurements for a is obtained. In this case, however, t can take 
any value. This trajectory is called “Pomeron trajectory”. The exchanged object carrying the quantum 
numbers of the vacuum (C = P = 1) is called “Pomeron”. 

The meson exchange leads to a simple power behaviour s“ _1 of total hadron-hadron scattering cross 
sections. In 1992 Donnachie and Landshoff 1299 1 made the observation that hadron-hadron total cross 
sections can be described by two powers, one a ss 1.1 corresponding to Pomeron exchange and one 
a = 0.5 corresponding to the exchange of mesons like p, w, fi and ai- Therefore at low-energies, where 
the Reggeon exchange dominates, a behaviour a ~ s -05 is expected, while at high energies u ~ s 01 
should prevail. 

Such a behaviour was indeed experimentally observed. Figure 22 shows the total inclusive proton- 
proton cross section as a function of the centre-of-mass energy fs for various experiments. In the 
range 2-20 GeV the proton-proton cross section decreases as measured by fixed-target experiments, while 
at higher energies it increases towards larger fs. The difference between proton-proton and proton- 
antiproton cross section decreases towards high energy; at about yfs ~ 30 GeV, the two cross sections are 
very similar. The rise of the proton-antiproton cross section for -y/5^30 GeV was first observed at the SPS, 
confirmed by Tevatron and LHC, where the ATLAS [300], CMS 13011 and TOTEM [302] experiments 
measured at fs = 7 TeV (TOTEM measured also at fs = 8 TeV |303j ). Data from cosmic-ray experi¬ 
ments 1304143061 are able to pro vide useful measurements at energies up to 60 TeV that are not covered 
by colliders. Overlayed in figure 22 as curve is a parameterisation nm based on the In 2 s behaviour that 
describes the data well at high energy. The power-law behaviour can also describe the present data 1308] . 
Also shown in figure 22 are recent measurements of the elastic proton-proton cross section [294] 13021 . 


9.4 The Underlying Event and Multi-Parton Interactions 

An important aspect of hadron-hadron collisions, neglected in our discussion so far, is the observation 
that additional mostly low-energetic (soft) particles are produced. The simple picture of a hard scatter¬ 
ing between one parton from each hadron with subsequent showering and hadronisation at least needs 
to be complemented by a treatment of the beam-beam remnants, which is beyond the capabilities of 
pQCD. In addition, the primary interaction may be accompanied by the production of further particles, 
a circumstance which is usually modelled in the form of multi-parton interactions (MPI). The concepts 
employed for MPI by Pythia and Herwig++, for example, are described in Refs. [3091431 ll and I312j . 
respectively. A very useful overview of general-purpose MC event generators including all aspects of event 
generation for the LHC is given in Ref. [42] . 

An unambiguous association of a specific particle to the reaction from which it originates is impossible. 
Therefore, the extra activity in a hadron-hadron collision cannot be uniquely separated from effects like 
initial-state or final-state radiation. One thus conventionally defines the “underlying event” (UE) as 
everything except for the hard process of interest. Two strategies are applied to measure the UE: The 
first evaluates events that are collected requiring only minimalistic trigger conditions such that even very 
soft collisions are registered. These so-called “minimum bias” (MB) measurements are supposed to be 
dominated by underlying-event physics. 

In the second approach, reactions are considered that have a high rate for the production of a leading 
particle or jet with reasonably high px (for triggering purposes). The event is then separated into hard and 
soft contributions through a physically motivated definition of phase-space regions that are dominated by 
either the hard or soft component of a collision. Traditionally, this is achieved by geometrically subdividing 
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Fig. 23. (a) Sketch of the geometrical subdivision of events in azimuthal angle for UE measurements, (b) Average 
density of charged particles in the transverse region as a function of the leading-jet or leading-track px, or the px 
of the reconstructed Z boson. (Adapted from Ref. IS 13V . ) 


the event into different regions (“towards”, “away”, and “transverse”) with respect to the direction of 
the triggering particle or jet (figure [23]( a)). Reactions that are typically chosen for this second approach 
are dijet or Drell-Yan lepton-pair production. Implicitly, it is assumed that the soft particle production 
is uncorrelated to details of the hard process—an assumption that can be tested experimentally. 

As an example for such a strategy figure [23}) b) shows a measurement by ATLAS (3131 of the charged- 
particle density in the transverse region as a function of the transverse momentum that is chosen as scale 
of the event hardness. Three event types are compared: minimum bias, dijet production, and inclusive Z 
production, with event scales chosen to be the px of the leading track, the leading jet, and the Z boson, 
respectively. An approximately continuous transition from MB to jet events is visible in the inset plot. 
In the comparison between jet and Z events of the soft-particle density in the transverse region, which is 
attributed to the UE, differences are exhibited, which are even more pronounced in the average summed 
Pt density of all tracks (not shown). These can be explained in terms of a bias between the definitions 
of the event scale. In the case of the jet events further jets with transverse momenta larger than the px 
of the leading jet are forbidden by definition, while for the inclusive Z events additional jet activity with 
Px > px is possible. Taking this bias into account, the conjecture that the UE activity is independent of 
the hard process is supported by the ATLAS measurement and by similar measurements from CMS 1314] . 

Another quantity that is not predictable by perturbative methods is the average density of charged 
particles produced in MB events for example at central rapidity. Corresponding measurements usually 
must be performed early on at the beginning of a data-taking period, when the instantaneous luminosity is 
small and the probability of multiple hadron-hadron collisions in the same or adjacent bunch crossings, the 
so-called pile-up, is negligible. It should be emphasised that such measurements are of utmost importance 
for a precise understanding of LHC data as they are used to model individual pile-up collisions. At later 
stages of the data-taking and with high instantaneous luminosity, the occurrence of up to 10-100 of such 
pile-up collisions simultaneously to the hard reaction of interest demands that these MB-like collisions 
are very well understood and modelled in simulations. Otherwise it becomes very difficult to precisely 
calibrate jet energies (see Chap. 3 of Ref. ®) or to estimate the background to searches for new physics. 
The charged-particle density is also of much interest for heavy-ion physics, and its dependence on the 
centre-of-mass energy is shown in Chap. 9, Fig. 9.1, of Ref. |Tj for various combinations of colliding 
hadron-hadron beams. 


9.5 Double-Parton Scattering 

The basis for all perturbative results in the previous sections (see in particular section is the assump¬ 
tion that collinear factorisation holds and that the cross section of interest can be calculated following 
equation Q. To some extent it is a great fortune that such a simple picture with only one interacting 
parton from each colliding hadron and only longitudinal degrees of freedom describes such a huge amount 
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of measurements. Partially this is due to the fact that one is mostly interested in inclusive processes at 
large px- Looking into the regime of small transverse momenta, however, it becomes apparent that this 
picture must fail, since otherwise the relevant cross sections would grow beyond all limits for px —> 0. 
They would thus violate unitarity, i.e. become even larger than the measured total cross sections. So, 
somehow, the cross sections have to be regularised. Realising that at px values of the order of a few 
GeV, protons are probed at fractional momenta of 10' 3 or below, one can interpret the consequently 
large parton densities in such a way that more than one parton per hadron participates in a scattering 
process. In terms of the MPI, as explained above, this concept has already been successfully applied to 
model the part of the soft-particle production that manifests itself in the form of the underlying event. 
The possibility of two hard or semi-hard reactions with two partons from each of the colliding hadrons 
participating in their own parton-parton scatter has not yet been considered, though. 

This type of process is called “double-parton scattering” (DPS), in contrast to the single parton 
scattering (SPS) discussed previously. Double-parton scattering is one possibility to dampen the rise 
of the low-px SPS cross sections. In a simplified form, the impact of double-parton scattering on the 
production rate for a final state A + B can be parametrised by one effective cross section that can be 
written as 


C^efF — 


ma A -a B 


2 


a 


dps ’ 

A+B 


(7) 


where a a and ctb are the cross sections for the independent processes A and B , <t^5b is the cross 
section for processes A and B to occur simultaneously, and m is a symmetry factor that is equal to unity 
for indistinguishable final states and two otherwise. In this simple picture, <r e ff is a measure of the overlap 
in transverse size of the parton distributions in the colliding hadrons. A naive geometrical interpretation 
leads to an estimate of a e g ss 7rR 3 ~ 50 mb, where R p is the proton radius. 

Measurements have been performed for a number of processes A and B: 4-jet [315U3161 . 7+3-jet 1317b 
El, and W+ 2-jet [320 1 32 ll production. In each case, the second independent scattering reaction is cho¬ 
sen to have a dijet final state, because dijet production occurs at by far the highest rates. To differentiate 
double-parton scattering from normal high-multiplicity single parton scattering events, the particular 
properties of DPS events are exploited: the px of the jet pair (both jet pairs in 4-jet production) should 
be balanced and the azimuthal angle between the vectors of the final-state objects of the leading scat¬ 
tering and the lower-px jet pair should be uncorrelated, i.e. randomly oriented. Specialising to VF+2-jet 
events, observables sensitive to DPS can then be defined as 


Ipt(ji) +Pt(J2)| 
Ipt(j'i)| + Ipt(j2)| ’ 


AS = arccos 


( PT(p,£r s )+p T (ji,j2) \ 

l v |p T (p,^ isS )| + |Px (jl j J2 ) | J 


where px are the transverse momentum vectors of the jets j\ , j 2 , or of the jet pair (jj, j 2 ) and the decay 
products of the W boson, one muon and missing transverse energy E™ lss representing the undetected 
neutrino. Figure [dlja) shows, as an example, the measurement of Z\ rel px by CMS [32lj together with the 
composition of single parton scattering background and double-parton scattering signal, which peaks at 
small px differences, as determined in a template fit. Figure [24] (b) summarises determinations of a e s from 
various experiments as a function of the centre-of-mass energy. The “Corrected CDF” point includes an 
imoroved estimate [3221!223j for the fact that eouation 0 implies processes A and B to be inclusive 
such that also higher terms than double scatterings are parametrised by a e g. The bias of exclusive event 
selections, when applying equation ([8]) in the W+ 2-jet case, 


_ +0-jet 

^eff — DPS * ^2-jet i 
a W + 2-jet 


( 8 ) 


has been considered in the presented analysis from CMS m and a similar measurement by ATLAS 13201 . 
Otherwise, a e g would not be independent any more of the choice of scattering processes A and B. 

The validity of the presented simplified theoretical approach anyway relies on a couple of additional 
assumptions. For a recent critical overview see Ref. [324]. 
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Fig. 24. (a) Template fit of SPS background and DPS signal to data for the zi rel pT distribution in W+ 2-jet 
events as measured by CMS. (b) Extracted effective cross section a e ff as a function of the centre-of-mass energy 
as measured by various experiments. (Adapted from Ref. \321H . ) 


10 Summary and Outlook 

The excellent performance of the LHC machine and of the four major detectors around the ring made it 
possible to carry out precise measurements at very high particle, i.e. jet transverse momenta, high invari¬ 
ant masses, and for very complex event topologies with large particle multiplicities. The high luminosities 
and the centre-of-mass energies of 7 and 8 TeV reached for proton-proton collisions are unprecedented. 
These conditions allowed to assess the validity of the QCD theory in the TeV regime and to stress-test 
our current understanding of particle-production processes at such high energies. Many measurements 
in particular of vector-boson and jet-production processes reached an experimental accuracy of a few 
percent uncertainty only—true QCD precision tests. 

To meet the experimental challenges, theoretical predictions are needed that match the precision 
and complexity requirements. Accordingly, enormous efforts have been undertaken to provide better and 
better theoretical predictions accompanied and driven by exciting new theoretical developments and 
insights. This includes largely automated NLO and many new NNLO calculations, improved analytical 
and parton-shower resummation techniques, the merging (i.e. matching) of parton showers and exact 
matrix elements, and attempts to gain a better understanding of the low-energy (i.e. non-perturbative) 
regime of the strong interaction. A major role is played by the Monte Carlo event generators that have 
been developed and optimised for the analysis and interpretation of LHC measurements. 

Calculations at NNLO in the strong coupling provide predictions for inclusive production cross sections 
with an estimated uncertainty, e.g. due to missing higher orders, of a few percent only. Such precision 
predictions were successfully tested against the LHC data. Often the dominant theoretical uncertainty 
originates from the knowledge of the proton PDFs. Also in this field LHC measurements have already 
contributed to a refined understanding of the proton structure. 

Around the year 2002, two theoretical break-throughs have revolutionised the field of Monte Carlo 
generators: first, the exact matching of 2 -4 2 production processes evaluated at NLO accuracy with 
parton showers, and second, the consistent merging of tree-level matrix elements with varying final-state 
parton multiplicity with parton showers. These two innovative concepts enable fully exclusive predictions 
at particle-level with an inclusive production rate accurate to NLO or the first few hardest emissions being 
modelled through exact real-emission matrix elements. Consequently, Monte-Carlo simulations based on 
matching or merging of matrix elements and parton showers have become the standard tools in any LHC 
analysis. Driven by the enormous progress in the evaluation of NLO cross sections (which are largely 
automated by now), there appeared many multi-leg NLO calculations that can directly be matched with 
parton showers. 

The currently emerging new standard are simulations based on exact NLO calculations with increasing 
parton multiplicity all matched with parton showers and combined into an inclusive description of the 
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considered production process. However, this is not the end of the development: first prototypes of exact 
NNLO calculations matched with parton showers have recently been presented. These exciting theoretical 
developments pave the way to a further scrutiny of our understanding of QCD in the upcoming high- 
energy and high-luminosity LHC runs. In consequence, they will significantly boost the sensitivity and 
prospects to find potential new physics at the LHC. 
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